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Chapter  1   INTRODUCTION 
1 .  1  Overview  of  MCDM 

Decision  making  is  a  complex  process  which  involves  multiple, 
conflicting,  and  noncommensurab  1  e  decision  criteria.  Problems 
encountered  in  public  and  private  sectors  are  always  given  in 
multiple  criteria  format.  For  example,  we  may  select  an 
automobile  that  consumes  less  fuel  and  costs  less  than  $10,000. 
But  often  we  worry  whether  the  automobile  is  the  safest,  the  most 
prestigestious .  the  most  comfortable,  and  the  highest  resale  value 
after  5  years.  City  planners  may  find  an  urban-development 
project  attractive  since  it  creates  lower  taxes.  However,  they 
may  be  concerned  that  if  the  project  creates  a  better  environment, 
more  public  services,  and  higher  job  opportunity.  A  business 
executive  may  favor  a  strategy  that  boosts  its  stock  price  and 
share  of  market,  but  he  may  worry  if  this  strategy  is  accompanied 
by  a  good  company  image,  good  labor  relationship,  and  favorable 
long  term  growth  potential.  Examples  like  these  can  go  on 
forever.  Wherever  they  appear.  all  the  examples  share  the  same 
characteristics:  multiple,  conflicting,  and  noncommensurab  1  e 
decision  criteria. 

1.2   Purpose  Of  This  Study 

The  techniques  the  DM  uses  to  reach  a  decision  are  intuitive 
judgement  and  mathematic  computation.  To  make  decisions  based  on 
a  human  being's  intuitive  judgement  are  the  most  widely  practiced 
decision  process,  especially  in  the   highest  management  level. 


However,  making  decisions  based  on  intuitive  judgement  can  be 
dangerous  and  unreliable  since  real  world  problems  are  usually  out 
of  any  individual's  comprehes i on.  The  second  technique,  on  the 
other  hand,  helps  the  DM  make  rational  decisions  by  a  systematic 
analysis  of  the  possible  courses  of  action.  In  view  of  the 
importance  of  the  second  technique,  this  study  will  be  focusing 
on  the  mathematic  formulation  of  MCDM  problems  and  their  solution 
methods  so  that  the  DM  can  utilize  this  study's  results  to  make 
rational  decisions. 


1.3   MADM  vs.  MODM 

The    problems    of    MCDM    are   widely    diverse.       However,    they   are 

broadly     classified     by     Hwang     and     Yoon      [1981]     among     other 

researchers,     into    two    categories:     Multiple    Attribute    Decision 

Making   (MADM).    and   Multiple   Objective   Decision    Making    (MODM)         The 

following    statement     from    Hwang    and    Yoon     [1981]     can    best    describe 

MODM    and   MADM: 

"...  Multiple  Objective  Decision  Making  (MODM)  is  not 
associated  with  the  problem  where  the  alternatives  are 
predetermined.  The  thrust  of  these  models  is  to  design  the 
"best"  alternative  by  considering  the  various  interactions 
within  the  design  constraints  which  best  satisfy  the  DM  by 
way  of  attaining  some  acceptable  levels  of  a  set  of  some 
quantifiable  objectives.  The  common  characteristics  of  MODM 
methods  are  that  they  possess:  (1)  a  set  of  quantifiable 
objectives;  (2)  a  set  of  well  defined  constraints.  (3)  a 
process  of  obtaining  some  tradeoff  information,  implicit  or 
explicit,  between  the  stated  quantifiable  objectives  and  also 
between  stated  or  unstated  nonquantifiable  objectives. 
The  distinguishing  feature  of  the  MADM  is  that  there  is 
usually  a  limited  (and  countably  small)  number  of 
predetermined  alternatives.  The  alternatives  have  associated 
with  them  a  level  of  the  achievement  of  the  attributes  (which 
may    not    necessarily    be    quantifiable)    based    on   which    the    final 


decision  is  to  be  made.  The  final  selection  of  the 
alternative  is  made  with  the  help  of  inter-  and  intra- 
attribute  comparisons.  The  comparisons  may  involve  explicit 
or  implicit  tradeoff...". 

Hwang  and  Yoon  also  provide  a  table  for  the  easy  comparison   of 

MADM  and  MODM  problems. 


Table  1.1    MADM  vs.  MODM 


MADM  MODM 

Criteria  (defined  by)   Attributes  objectives 

Objective  Implicit  Explicit 

(ill  defined) 

Attribute  Explicit  Implicit 

Constraint  Inactive  Active 

(Incorporated  into 
attributes  ) 

Alternative  Finite  number.        Infinite  number. 

di sere t e ( pr es cr ibed)    Continuous 

(emerging  as  process 
goes  ) 

Interaction  with  DM     Not  much  Mostly 

Usage  Selection/Evaluation     Design 


To  cope  with  the  diversified  MCDM  problems.  researchers 
developed  solution  methods  for  the  two  categories.  MADM  and  MODM. 
respectively . 

1.3.1   MADM  Solution  Methods. 

The  earlier  works  on  MADM  are  from  many  different  disciplines: 
management  science,  economics,  psychome t r i cs .  marketing  research, 
applied  statistics,  decision  theory,  and  so  on  [Hwang  and  Masud. 
1981],  Each  discipline  confronted  multiple  criteria  in  quite 
different  situations.  Consequently  each  area  had  developed 
method(s)  for  its  own  particular  usage.  Hwang  and  Yoon  [1981] 
collected  and  classified  17  methods  based  on  the  information 
requirements  from  a  DM.  Figure  1.1  depicts  that  classification 
scheme  in  great  detail.  Each  method  involves  its  own  types  of 
underlying  assumptions,  information  requirement,  and  evaluation 
principles.  Some  methods  may  not  fully  satisfy  the  underlying 
assumptions  (e.g.  ,  complete  independence  among  attributes  for 
simple  additive  weighting  method),  and  some  may  have  more  tedious 
procedures  (e.g.  ,  constructing  indifference  curves  for 
hierarchical  tradeoffs)  than  others.  There  are  also  methods 
fitted  for  a  specific  situation  (e.g.  ,  maximin  for  the 
pessimistic  decision  situation.  maximax  for  the  optimistic, 
disjunctive  method  for  the  selection  of  a  specialist).  Even  with 
the  same  preference  information,  each  method  evaluates  the  problem 
from  a  different  point  of  view  (e.g.  .  with  a  set  of  weights, 
weighted  average  outcomes  for  simple  additive  weighting,  relative 
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closeness  to  the  ideal  solution  for  TOPSIS.  satisfaction  to  a 
given  concordance  measure  for  linear  assignment  method  and 
ELECTRE) . 

It  is  argued  by  Zeleny  [1982]  that  alternatives  that  are 
closer  to  the  ideal  are  preferred  to  those  that  are  farther.  To 
be  as  close  to  the  perceived  ideal  is  the  rationale  of  human 
choice.  From  this  argument,  we  found  that  TOPSIS  which  was 
developed  by  Hwang  and  Yoon  [1981]  may  be  the  the  best  among  all 
the  MADM  methods.  The  reasons  are  that  it  considers  not  only  the 
positive  ideal  point  but  also  the  negative  ideal  point,  it  is  easy 
to  calculate  and  to  program  (as  a  matter  of  fact.  TOPSIS  Computer 
software  has  been  developed  at  Kansas  State  University),  and  its 
power  to  handle  large  problems.  In  this  study.  TOPSIS  will  be 
further  improved  to  incorperate  various  normalization  equations 
and  distance  measurements. 


1.3.2    MODM  Solution  Methods. 

Like  MADM  solution  methods,  the  works  on  MODM  are  from 
various  research  areas.  Each  method  has  different  types  of 
assumptions,  different  information  requirements,  and  are  even 
designed  for  a  particular  type  of  problem.  Hwang  and  Masud  [1979] 
had  collected  and  classified  19  methods.  Figure  1.2  depicts  that 
scheme  of  classification  in  detail. 

The  first  class  requires  no  information  from  the  DM  once  the 
problem  constraints  and  objectives  have  been  defined.  The  analyst 
makes  assumptions  about  the  DM's  preference  and  presents 
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solution(s)  to  the  DM.  The  second  class  assumes  the  DM  has. 
consciously  or  subconte lous  1  y .  a  set  of  goals  to  achieve  and  the 
goals  set  is  given  to  the  analyst  before  the  formulation  of 
mathematic  model.  The  third  class,  also  known  as  interactive 
methods,  requires  greater  DM  involvement  in  the  solution  process. 
The  interactive  takes  place  through  a  DM-analyst  or  DM-computer 
dialogue  at  each  iteration.  The  tradeoff  or  preference 
information  given  by  the  DM  after  each  iteration  is  used  for 
determining  a  new  solution.  The  DM.  hence.  actually  gains 
insights  into  the  problem.  Finally,  the  fourth  class  does  just 
one  thing  -  to  determine  a  subset  of  the  complete  set  of 
nondominated  solutions  to  the  MODM  problem.  In  doing  so.  it  deals 
strictly  with  the  physical  constraints  and  makes  no  attempt  to 
consider  the  preference  of  the  DM.  The  desired  outcome,  however. 
is  to  narrow  down  the  possible  courses  of  actions  and  makes  the 
selection    of    the    preferred    course    of    action    easier. 

None  of  the  processes  are  superior  to  the  others,  because 
one  method's  advantages  is  another's  shortcomings  and  vice  versa. 
Therefore,  two  algorithms  were  proposed  to  tackle  MODM  problems. 
The  first  method  has  less  DM  involvement  and  is  called  TOPSIS- 
MODM.  The  second  one  requires  great  involvement  of  DM  and  is 
called    Interactive    Sequential    Goal    Programming    II    (ISGP    II). 


1.4         Theses    Structure 

Three   algorithms    for   complex   MCDM   problem  were   attempted   in 
this    study.       Chapter     1     briefly    reviews    the    MCDM    and    its    solution 


methods.  Chapter  2  describes  a  MADM  algorithm  called  TOPSIS  II 
which  is  a  modified  version  of  TOPSIS  developed  by  Hwang  and  Yoon 
in  1981.  Chapter  3  presents  a  MODM  algorithm  called  TOPSIS-MODM 
which  utilizes  the  structure  of  the  global  criteria  method  and  the 
concept  of  TOPSIS-MADM.  Chapter  4  presents  an  interactive  MODM 
algorithm  called  ISGP  II  which  is  a  refinement  of  Hwang  and  Yoon's 
ISGP  [1982].   Finally  the  conclusions  are  given  in  Chapter  5. 


Chapter    2         TECHNIQUE    FOR    ORDER    PREFERENCE    BY    SIMILARITY    TO    IDEAL 

SOLUTION    II     (TOPSIS    II) 
2  .  1       Introduction 

Hwang  and  Yoon  [1981]  developed  the  Technique  for  Order 
Preference  by  Similarity  to  Ideal  Solution  (TOPSIS)  based  upon  the 
concept  that  the  chosen  alternative  should  have  the  shortest 
distance  from  the  ideal  soultion.  A  .  and  the  farthest  from  the 
negative-ideal  solution.  A".  That  is.  for  a  MADM  problem  with 
each  attribute  having  a  monotonical 1 y  increasing  (or  decreasing) 
utility,  a  preferred  alternative  should  be  closest  to  the  ideal 
solution  --  which  is  composed  of  all  the  best  attribute  values  -- 
and  farthest  from  the  negative-ideal  solution  --  which  is  composed 
of  all  the  worst  attribute  values.  However,  it  is  not  easy  to 
find  a  preferred  alternative  which  satisfies  the  foregoing 
situation.  For  example,  in  Fig.  2.1.  an  alternative  Aj  has  shorter 
distances  to  both  the  ideal  solution  and  the  negative-ideal 
solution  than  alternative  A2  does.  It  is  then  very  difficult  to 
justify  the  selection  of  A  j .  To  cope  with  this  difficulty. 
TOPSIS  measures  relative  closeness  to  the  ideal  solution  rather 
than  absolute  distances  to  both  ideal  and  negative-ideal  solutions 
to    determine    the    preferred    alternative. 

It  was  found  later  by  Yoon  [1986]  that  even  with  relative 
measurements  of  distance,  conflicting  preference  alternatives  may 
actually  arise  if  different  parameter  values  are  used  in  the 
distance  function  and/  or  normalization  function  of  TOPSIS.  To 
reconcile  that  conflict,  modification  of  TOPSIS  is  made  by 
introducing    the    credibility    of    a    distance    function,     d    . 
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Figure  2.1  Distances  to  The  Ideal  and  Negative  Ideal 
Solutions  in  Two  Dimension  Space. 
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Normalisation  of  Attribute  Ratings^. 

For  multiple  attribute  decision  making  (MADM)  problems,  the 
measurement  units  that  each  attribute  possesses  are  different. 
This  makes  the  direct  comparisons  among  attributes  impossible. 
Hence,  the  normalization  of  the  attribute  ratings  is  a  necessary 
step  in  solving  a  MADM  problem. 

Some  different  normalization  equations  have  been  developed  by 
Zeleny  [1982],  Hwang  and  Yoon  [1981],  and  Nijkamp  and  Van  Delft 
[1977].  The  combination  of  these  equations  is  done  by  Yoon  [1986] 
using  the  L  metric.  If  one  defines  the  normalized  scale  as  the 
ratio  between  individual  and  combined  distance  from  the  origin, 
the  comparable  scale  of  alternative  i.  Aj  according  to  attribute 
j.  Xj  is  defined  as: 


(q)  -  X 


1J  I    <  ?   Xi 


q  >i/q 


(2.1) 


where  q  can  be  any  integer  between  1  and  oo.   For  q  =  oo.  each 
attribute  is  compared  by  its  largest  attribute,  i.e.: 


rjj  (oo)  =  Xjj  /  max  {  Xjj  > 


(2.2) 


Three  parameters,  q  =  1.  2.  and  oo  are  widely  used.  Expressions 
(2.1)  and  (2.2)  will  be  used  to  transform  the  various  attribute 
dimensions  into  nond imens iona  1  attributes,  allowing  the  direct 
comparisons  across  the  attributes. 
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Measures  of  Distance^. 

The  notion  "closeness"  can  be  interpreted  geometrically  as  the 
distance  between  two  points  in  a  n-dimension  space,  e.g.  given 

o2 xn2).  the 


points  x1  ^Ix^.xj1 x/  )  and  x2  -  (x^ 

distance  d  between  them  is  found  to  be: 


d„  =  {  Z   W, 


V-  M2,P  > 


1/p 


(2.3) 


where  N.  are  weights;  1  <  p  <  oo. 

For  p  =  oo.  expression  (2.3)  is  rewritten  as 


«ax  {  Wj  (  x.- 


Xj2)   > 


(2.4) 


We  can  also  see  why  the  values  p  -  1.  2.  and  oo  are 
especially  operationally  important:  dj  (  the  Manhattan  distance  or 
city  block  distance  )  and  d2  (the  Euclidean  distance)  are  the 
longest  and  shortest  distances  in  the  geometric  sense:  dQ0  (the 
Tchebycheff  distance)  is  the  shortest  distance  in  the  numerical 
sense . 

Expressions     (2.5)    and    (2.6)    are    used    to    determine    the    separation 
measure    between    an    alternative    i    to    the    ideal    solution. 


dp*(i)    -    (    ?       W/    |r/    -     riji    P    )X/P 
for    1    <    p    <    oo       ,    and 

d    *(i)     ■     (     max    Wj     (r j"    -    r j j )     }  ,       P    -    oo 


(2.5) 


(2.6) 
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Similarly,  the  separation  measure  between  alternative  i  to  the 
negative  ideal  solution  is  written  as: 


dp"(i)  =  {  Z      W 


rij  I"  > 


P  i  1/P 


r  j 

d0o"(il  =  {  max  Wj  I  rj"  -  i-jj  I  ) 


(2.7) 


(2.8) 


Conflict  arises  when  selection  of  the  best  alternative  based  on 
d  'is  different  from  that  of  dp~.  To  resolve  t h i sco n f 1 i c t .  the 
relative    closeness    of    Aj    with    respect    to    A      is    introduced    and    is 


written    as  : 


D„(i)    =    d    ~U)    /    {    d    *(i)    ♦    d    "(1)    } 


(2.9) 


where    i    =    1.2 m:     p    ■    1.2.     and    oo 

* 

It    is    clear    that    0    <    D(i)    <     1    where    D(i)    -    1     if    Aj    =    A    ,     and 
D(i)   ■  0   if  Aj   ■  A    .      An  alternative  Aj    is   relatively   closer   to  A 
as    D(i)    approaches     1. 

The  concept  of  relative  closeness  to  A  resolves  the  selection 
problem  of  the  reference  point  between  A  and  A".  TOPSIS  still 
renders  up  to  nine  different  solutions  or  nine  sets  of  alternative 
rankings  by  changing  the  p  and  q  parameters  in  both  the 
normalization  and  separation  computation  steps.  A  crude  resolution 
among  conflicting  rankings  can  be  made  by  introducing  the  relative 
reliability    rate    of    a    distance    (dp)    measurement. 
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&£££E5!iaSli2!l  of  Cl£dibil.  j.  ty_  of  A  Distance  Function^.  dp: 

The  physical  property  of  d  plays  a  key  role  in  determining 
the  credibility  of  itself.  That  property  can  be  best  described 
using  indifference  curves.  An  indifference  curve  is  composed  of 
many  points.  Any  point  on  the  same  indifference  curve  has  equal 
distance  from  the  given  origin.  The  longer  the  curve,  the  more 
points  it  can  contain.  That  is.  more  points  can  be  expressed  by 
the  same  distance  with  a  longer  indifference  curve.  For  example. 
in  Fig.  2.2,  there  are  three  loci  on  a  two-dimensional  space  with 
p  -  1,  2.  and  oo.  Each  locus  stands  for  an  indifference  curve. 
Notice  that  the  locus  appears  farther  from  the  origin  as  p 
increases.  This  shows  that  L  becomes  longer  as  p  increases. 
Consequently,  more  points  can  be  contained  in  the  locus  with  a 
larger  p  or  more  points  can  be  expressed  by  the  same  dp  with  a 
larger  p.  We  conclude  that  the  distance  function  d_  becomes  less 
specific  (or  less  reliable)  as  p  increases. 

If  we  view  the  distance  function  as  a  measurement  system. 
such  as  a  grading  system,  temperature  reading  scale,  and  so  on. 
the  foregoing  argument  will  be  more  clear.  A  measurement  system 
is  more  credible  when  it  has  more  classes.  In  another  words,  a 
more  credible  measuring  system  has  a  smaller  class  size.  For 
example,  the  Fahrenheit  scale  is  more  credible  than  the  centigrade 
scale  because  the  former  divides  the  temperature  range  between 
freezing  and  boiling  into  180  degrees  (classes)  while  the  latter 
divides  this  range  into  only  100  degrees  (classes).   A  temperature 
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Figure  2.2  Indifference  Curves  with  Different  d 
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expressed  by  20°C  can  be  denoted  by  67°F.  68°F,  or  69°F  which  is 
more  precise  and  reliable.  Since  a  class  contains  only 
homogeneous  elements  with  respect  to  its  concern,  an  indifference 
curve  is  a  class  as  far  as  the  measurement  of  distance  is 
concerned.  Hence  the  length  of  an  indifference  curve  can  be 
viewed  as  a  class  size.  Recall  that  the  length  of  an  indifference 
curve  (class  size)  increases  as  p  increases.  Hence,  we  conclude 
the   credibility  of  d   decreases  as  parameter  p  increases. 

Since  the  relative  reliability  of  d  is  proportional  to  the 
average  of  dD,  E(d_).  the  quantification  of  the  credibility  of  dp 
is  done  by  using  E(d_).  The  procedure  to  obtain  the  values  of 
E(d_)  are  as  follows. 

1.  Generate  (  n  *  30000)  random  numbers  whose  range  is  (  0.0. 
1.0)  where  n  is  the  dimension  (e.g.  the  number  of  attributes  or 
the  number  of  alternatives)  of  the  problem. 

l.a.  Take  the  summation  of  these  random  numbers  and  divide 
it  by  30000  to  get  E(dj). 

l.b.  Take  the  square  of  each  random  number  and  add  all  the 
results  together.  After  the  summation,  take  the  square-root  of 
that  summation,  divide  it  by  30000.  and  we  get  E(d2). 

I.e.  Group  the  random  numbers  into  30000  n-element  pairs. 
For  example,  for  a  2-dimensional  problem,  the  first  random  number 
goes  with  the  second  random  number,  the  third  random  number  goes 
with  the  fourth  random  number,  and  so  on.  Compare  within  each 
group  and  find  the  largest  number  for  that  group,  and  add  those 
largest  numbers  together.   Take  that  summation  and  divide   it  by 
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30000:   we  get  E(d00).   Table  2.1  displays  the  results  of  the 
foregoing  procedure. 

With  all  the  E(d  )  values,  the  relative  credibility  of  each  dp 
under  the  three  distance  system  is  written  as: 


Cp=    E(dp)     /     {     Etdj)     +    E(d2)     +    E(d00)     >.     p=    1.     2.     oo. 

(2.10) 
Table     2.2    displays     the     results     for     all     the     Cp     values. 
(Appendix    A    shows    a    BASIC    program    for    generating    E(dp)    and    Cp). 


E.£°a£iIiS£  So.lutionsj_ 

As  expected.  C,  increases  very  rapidly  when  the  number  of 
dimensions  increase.  From  a  purely  mathematlc  view,  the  use  of  Cj 
in  determining  the  final  solution  is  trustworthy.  However,  if  we 
consider  three  distance  measurements  simultaneously,  the  unified 
function    is    given    as: 

(2.11) 


Notice  that  C  values  are  the  weights  In  this  linear  combination 
form  and  dj  is  the  most  heavily  weighted  distance  measurement 
among  all  dp. 

Similarly,  normalized  attribute  ratings  get  less  specific  when 
q  increases.  Therefore,  the  same  set  of  C  can  also  be  applied  in 
the  unified  normalization  equation. 


■ti    =  Cj  rjjll)  +  C2  eii{Z)    *   CQO  rij(oo). 

/ 


(2.12) 
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Table  2.1  Mean  Values  of  Distance  Measurements 


Dimension 


E(d,  ) 


Eld,) 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

la 

20 


1 .0026 
1 .5046 
1 .9954 
2.4913 
2 .9956 
3.4975 
4 .0009 
4 .5043 
4.9979 
5.5047 
6.0010 
6.4949 
6 .9964 
7 .4999 
8.0064 
8 .4895 
9.0040 
9.5045 
9 . 9857 


0 . 7671 
0.9635 
1 . 1199 
1 . 2582 
1 .3873 
1 .5036 
1.6115 
1 . 7130 
1 .8063 
1 .8984 
1 .9820 
2.0641 
2. 1434 
2 . 2201 
2 .2954 
2.3638 
2.4364 
2.5031 
2 .5662 


0.6681 
0.7520 
0. 7990 
0.8322 
0.8571 
0.8756 
0.8885 
0.9006 
0 .9093 
0.9180 
0.9228 
0 .9290 
0.9329 
0 .9374 
0 . 9415 
0.9443 
0 . 9477 
0.9497 
0. 9522 


49 
5  0 


24 . 5066 
25 .0146 


4 . 0338 
4 .0765 


0 . 9800 
0. 9804 
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Table  2.2  The  Relative  Credibility  of  Each  Distance  Function 


Dimension 
1 
2 
3 
4 
5 


1 

.3333 

.4113 
.4673 
.5098 
.  5437 
.5717 
.5951 
.6154 
.6328 
.6479 
.6615 
.  6738 
.  6845 
.  6946 
.  7037 
.7121 
.  7196 
.  7268 
.7335 
.7395 


.3333 
.4137 
.2992 
.2861 
.2747 
.2647 
.2559 
.2479 
.2407 
.2342 
.2281 
.2225 
.2175 
.2128 
.2083 
.  2042 
.2004 
.  1967 
.  1932 
.  1900 


oo 
.3333 

.  2741 

.  2335 

.2041 

.  1816 

.1636 

.  1490 

.1367 

.  1265 

.  1179 

.1103 

.  1036 

.0979 

.  0926 

.  0880 

.0837 

.  0800 

.0765 

.0733 

.0705 


9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 


49 
50 


.8302 
.  8318 


.  1366 
.  1356 


.  0332 
.0326 
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The  only  difference  between  expressions  (2.11)  and  (2.12)  is 
that  the  criteria  used  in  determing  C  is  different.  For  (2.11).  C 
is  decided  by  the  number  of  attributes  while  C  in  (2.12)  is 
determined  by  the  number  of  alternatives. 
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2.2    Problem  Formulation  and  Algorithm 


Problem : 


D  = 


<21 
lll 

Kml 


X 


2 

x12 
x22 

xi2 

xm2 


Aj 

XU 

x2j 

Xij 

xmj 


Xn 

xln 

x2n 

xin 


AlsoiiiSisi 

ilSE  i-i.    Normaiization  of  Problem  Matrix  X 

Wj  [  tijl  q  )]  -  Wj  [  Xjj  /  {  Z   Xjj"  }1/q  ] 
where  q  can  be  any  value  between  1  and  oo  .  I.e..  1  <  q  <  oo. 
However,  the  most  commonly  used  q  values  are  1.2.  and  oo.   When  q 
=  oo.   each  attribute  is  compared  by  only  its  largest  attribute. 
That  is 

Wj  [  rijloo)]  =  Wj  [  Xjj  /  {  I   Xij00  >  1/0°   ] 

can  be  rewritten  as: 

Wj  [  rtj(oo)]  =  Wj  [  X^  /  max  {  Xtj  >  ] 

Ste2   2^_        Unification    °L   lllESS   S2E!SSlil£^    Matrices 

Different    q    values    result    i nd i f f e r en t nor ma  1 i z e d    matrices. 
The    combination    of    the    three    matrices    follows    the    formula: 
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ij  =  Cj  ri;i(n  ♦  C2  ri}{2)   *   C00  rjjlool 


where  Ci,  C2  and  C00  are  determined  by  the  number  of  alternatives. 
See  Table  3.2  for  C_  values. 


Steg  3j.   DiSElay.  The  Unified  Matrix 


Rij  *   A2 


'1  1 


R,-2 


sil     Ki 


*ml     nm2 


<li 
*2J 


<ij 


RmJ 


In 


Step.  4_^    Determination  of  Positive-ideai  Soiution  and 
Negative-ideai  Soiution 

R   «  {(max  Rj.j  |  j  E  Jl.lnin  Rlf  |  j  E  J')  |  i  -  1.2 m> 

•    <  Rx  •  R2  Rj  Rn   } 

R~  =  {(min  R ±  1  |  j  E  J). (max  Rjj  I  j  E  J')  |  i  -  1.2 m) 


-    f  Ra*.  R2 Rj 


•  Rn~  > 


where  J  is  the  higher  the  better  attribute,  and  J1  is  the   smaller 
the  better  attribute. 
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Step.  Ll  £al£!llat   ion  ofd   (i)  and  dp~(i) 


dp*(i)  -  {  Z   W/  I  Rj*  -  Rjj  lp  }  1/p 


I* 


d  "(i)  =  {  Z      W,P  I  RT  -  Rn  lp  >  1/p 

P         3-1    J      J 

where  p  can  be  any  integer  between  1  and  oo,  i.e.  1  <  p  <  oo. 
However,  the  most  commonly  used  p  values  are  1.2.  and  oo  .  When  p 
=  oo.  the  formulations  of  d  *(i)  and  dp~(i)  can  be  rewritten  as 


dp  (i)  ■  max  Wf  I  Rj   -  Rjj  I 


d   (i)  =  max  Hj  |  Rj   -  Rj  i  I 


St eg  6^    Unification  of  d   and  d~ 

Different  p  values  result  in  different  distance  measurements. 
The  combination  of  d    can  be  done  by  using: 

d*(i)  -  Cx  da*(i)  +  C2  d2*fi)  ♦  CQ0  dQO*(i) 

Similarly,  the  combination  of  d  ~  can  be  done  by  using: 

d"(i)  =  Cj  dj"{l)  +  C2  d2"(i)  <•  C00  d00"(i) 

where     C j  .      C2.      and     C00     are     determined     by     the     number     of 
attributes.        See    Table    2.2    for    C       values. 

S£ep_    7^         Ranking   The    Preference    Order 

There    are    d      and   d"    values    for    each   alternative    (Aj).       If   d 
is    used,     one    chooses    the    alternative    with    the    smallest    d       as    the 
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best  alternative.   If  d   is  used,  the  alternative  with  largest  d 
is  thought  to  be  the  best. 

However,  conflict  occurs  when  the  ranking  orders  for  d   and  d 

are  different.   To  resolve  the  conflict  ,  a  unified  D  =  d"  /  {  d~ 

+  d  }  is  introduced.  The  alternative  with  the  larger  D  value  is 
thought  to  be  superior  than  the  one  with  the  smaller  D  value. 
Hence,  a  unique  ranking  order  is  obtained. 
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2.3  Numerical  Examples 

2.3.1  Selection  of  Fighter  Aircraft 

A  country  decided  to  purchase  a  fleet  of  jet  fighters  from  the 
U.S.  Pentagon  officials  offered  the  characteristic  information  of 
four  models  which  may  be  sold  to  that  country.  The  Air  Force 
analyst  team  of  that  country  agreed  that  six  characteristics 
(attributes)  should  be  considered.  They  are:  maximum  speed  (X.). 
ferry  range  (X2).  maximum  payload  (Xg),  purchasing  cost  (X4). 
reliability  (X5).  and  maneuverability  (Xg).  The  value  of  the  six 
attributes  for  each  model  (alternative)  are  given  in  Table  2.3. 

Table  2.3   Decision  Matrix  for  Aircraft  Selection: 


D  - 


l4 


Xl  X2  X3  X4  X5 

2.0  1500  20000  5.5  average 

2.5  2700  18000  6.5  low 

1.8  2000  21000  4.5  high 

2.2  1800  20000  5.0  average 


A6 
very  high 

average 

high 

average 


Before  the  calculation,  we  quantified  the  qualitative  data  using 
the  bipolar  scale  developed  by  MacCrimmon  [1968].  This  scaling 
system  is  shown  in  figure  2.3. 
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For  cost  attributes 


For  benefit  attributes 


very  high    1.0 

high        3.0 

average 

low 
very  low 


5  .  0 

7.0 

9.0 

10  .0 


0 

1 . 0  very  low 

3.0  low 

5 . 0  average 

7.0  high 

9.0  very  high 

10.0 


Fig.  2.3   Scaling  Systen 


The  decision  matrix  thus  becomes: 


D  « 


xl 

X2 

X3 

X4 

X5 

X 

2.0 

1500 

20000 

5.5 

5 

9 

2  .  5 

2700 

18000 

6.5 

3 

5 

1  .8 

2000 

21000 

4.5 

7 

7 

2  .  2 

1800 

20000 

5  .0 

5 

5 

Note  that  all  attributes  except  X4  are  benefit  criteria. 

We  assume  the  DM  assigns  the   relative   weight  to   each   attribute 

as:   W  «  (  0.30.  0.20,  0.10.  0.10.  0.10.  0.20  ) 


Stej)    1-        S°£!!!Siization    of    Deci.si.on    Matrix 

For     the     q     =     1     case,      expression     (2.1)     is     used     and     the 
normalized    weighted    decision    matrix    is    written    as: 


D    = 


0706 

.0375 

.  0253 

.0256 

.  0250 

.  0692 

0882 

.0675 

.0228 

.0302 

.0150 

.0384 

0635 

.  0500 

.  0266 

.  0209 

.0350 

.0538 

0776 

.0450 

.0253 

.  0233 

.  0250 

.0384 
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Similarly,  we  obtain  the  second  normalized  weighted  decision 
matrix  for  q  »  2 . 


1402 

.  0732 

.0506 

.  0506 

.0481 

.  1342 

1752 

.  1318 

.0455 

.0598 

.0298 

.0745 

1261 

.  0976 

.0531 

.0414 

.0674 

.1043 

1542 

.0878 

.  0506 

.0460 

.  0481 

.0745 

For  q  =  oo  .  expression  (2.2)  is  used  and  the  matrix  is  written  as: 


.2400 

.  1111 

.0952 

.0846 

.0714 

.0714 

.3000 

.2000 

.0857 

.  1000 

.  0429 

.1111 

.  2160 

.  1481 

.  1000 

.  0692 

.  1000 

.1556 

.2640 

.  1333 

.  0952 

.0769 

.0714 

.1111 

Step_  2^    Uailicaiiaa  °t  l!l£££  Normalized  Weighted 
5££ilion  Matrices 
The  combination  of  the  three  matrices  is  done  using  expression 
(2.12).  where  the  relative  reliability  coefficients.  C  ,  for  the 
four  alternatives  are  .5098.  .2861.  and  .2041  for  q  =  1,  2.  and  oo 
respectively.   The  unified  decision  matrix  is  written  as: 


246 

.0626" 

.0466 

.0447 

.04  10 

.1142 

.  1558* 

. 1126* 

.0420" 

.0527" 

.0246" 

.0634 

.  1122" 

.  0834 

. 0490* 

. 0365* 

.0574* 

.0888 

.  1371 

.0751 

.0466 

.0406 

.  0410 

.  0634 

28 


Ste£    3^         Determination    of    Pos i t ive-Idea  1     and    Negative-ideal 
Solutions 

The    positive-ideal  solution    is: 

R*    -     (     .1558.  .1126.        .0490.        .0365.        .0574.        .1142    ) 

The    negative-ideal  solution    is: 

R"    ■     (     .1122.  .0626.        .0420.        .0527,         .0246.         .0634     ) 

Step.  4^_    Caicuiation  Of  Separation  Measures 

For  the  p  =  1  case,  expressions  (2.5)  and  (2.7)  are  used  to 
calculate  d  (i)  and  d  ~(i)  respectively,  and  the  result  is  given 
in  Tabl e  2.4. 

Table  2.4  Separation  Measures  When  p  ■  1 

i        dj'd)  dj"(i) 

1  .1082  .0922 

2  .1068  .0936 

3  .0982  .1022 

4  .1299  .0705 

Similarly,  separation  measures  for  the  p  =  2  case  are  obtained 
and  are  shown  in  Table  2.5. 
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Table  2.5  Separation  Measures  When  p  =  2 


1 

.0617 

2 

.0630 

3 

.0583 

4 

.0680 

.0556 
.0663 
.0453 
.  0348 


For  the  p  =  oo  case,  expressions  (2.6)  and  (2.8)  are  used  to 
calculate  d  (i)  and  d  ~(i)  respectively,  and  the  result  is  shown 
in  Tabl e  2.6. 


Table  2.6  Separation  Measures  When  p  =  oo 
i  d   *<i)  d   "(i) 


.0500 
.0508 
.0436 
.0508 


.0508 
.0500 
.  0328 
.0249 


£££E  5-   Unification  Of  1".°  iSESJiation  Measures^  d    and  d_~ 

Unification  of  d  is  obtained  using  (2.11).  Similarly, 
unification  of  d~  is  obtained  using  (2.11).  For  six 
attributes,  C 1  -  .5717  .  C2  -  .2647, and  CQ0  =  .1636  are  used. 
The  results  are  given  in  Table  2.7 
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Table  2.7  Unified  d    &  d    and  Their  Rankings 


d  ( i )    ranks 


d  ( i )   ranks 


0863 

3 

0860 

2 

0787 

1 

1005 

4 

0757 

2 

0793 

1 

0769 

3 

0536 

4 

§l£E  §.+.   Banking  The  Preference  Order 

Notice  the  ranking  order  for  d  is  different  from  the  one  for 
d".  This  conflict  is  resolved  using  expression  (2.9).  The  values 
for  all  D's  and  the  final  ranking  order  is  shown  in  Table  2.8. 

Table  2.8  Final  Solution 
1       D(i)  Rank 


1 

.4673 

2 

.4797 

3 

.4942 

4 

.3478 

The  solution  is: 
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2.3.2   Helicopter  Evaluation  Problem 
A  company  currently  evaluates  five  helicopters  for  its  president. 
The  quantifed  data  matrix  is  given  as: 


xl 

X2 

X3 

X4 

X5 

X 

• 

145 

580 

1  .9 

2625 

3.5 

9 

190 

500 

2.  2 

2700 

3.0 

5 

150 

450 

1  .8 

2550 

2  .5 

9 

135 

620 

2.5 

2700 

4  .  5 

3 

170 

430 

2.0 

2600 

4  .0 

7 

All  the  attributes  except  Xg  are  benefit  criteria. 

We  also  assume  the  DM  assigns  the  relative  weight  to  each 

attribute  as  W  ■  (  1.0.  1.0.  1.0.  1.0.  1.0.  1.0.  1.0.  1.0  ) 

Step.  1^_    Normalization  of  Decision  Matrix 

For  the  q  =  1  case,   expression  (2.1)  is  used  and  the 
normalized  weighted  decision  matrix  is  written  as: 

.1835   .2248   .1827   .1992   .2000   .2727   .2414   .2800^ 


.2405   .1938   .2153 
.1899   .1744   .1731 


.2049   .1714   .1515   .1034   .2800 
.1935   .1429   .2727   .3103   .2000 


.1709   .2403   .2404   .2049 


.  2571 


.  0909 


.  1724 


.  2000 


.2152   .1667   .1923   .1973   .2286   .2121    .1724   .0400 

Similarly,  we  obtain  the  second  normalized  weighted  decision 
matrix  for  q  =  2. 
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D  = 


.4073 

.4977 

.4056 

.4454 

.4384 

.  5750 

.  5092 

.5735 

.5337 

.  4290 

.4697 

.4581 

.3757 

.3194 

.2182 

.5735 

.4213 

.3861 

.  3843 

.4327 

.3131 

.5750 

.6547 

.4096 

.3792 

.5320 

.5337 

.4581 

.5636 

.  1917 

.  3637 

.4096 

.4775 

.3690 

.4270 

.4412 

.  5010 

.4472 

.  3637 

.0819 

For  q  ■  oo,  expression  (2.2)  is  used  and  the  matrix  is  written  as 


D  - 


9355  .7600  .9722   .7778  1.0000  .7778  1.0000 

1.0000   .8065  .8800  1.0000   .6667  .5556  .3333  1.0000 

.7895   .7258  .7200  .9444   .5556  1.0000  1.0000  .7143 

.7105  1.0000  1.0000  1.0000  1.0000  .3333  .5556  .7143 

.8947   .6953  .8000  .9630   .8889  .7778  .5556  .1429 


S  t  efi  2j_     Unification  of  Three  Norma  1.  j.  zed  Weighted  Decision 
Matrices 
The  combination  of  the  three  matrices  is  done  using  expression 
(2.12).  where  the  relative  reliability  coefficients.  C  ,  for  the 
five  alternatives  are  .5437.  .2747.  and  .1816  for  q  =  1.  2,  and  oo 
respectively.   The  unified  decision  matrix  is  written  as: 


.3464  .4242  .3450  .4028  .3664  .4825  .4078  .4860 

.4540*  .3656  .4015  .4143  .3140  .2681  .1747"  .4860* 

.3584  .3291  .3268"  .3913"  .2817*  .4825  .5243*  .3472 

.3226"  .4534*  .4539*  .4143*  .4710"  .1608"  .2913  .3472 

.4062  .3154"  .3631  .3990  .4187  .3753  .2913  .0694" 
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Step_  3^   Determination  of  Positive-ideai  and  Negative^ideai 
Solutions 

The  positive-ideal  alternative  is: 
R*  -  (  .4540.  .4534.  .4539,  .4143.  .2817.  .4825.  .5243.  .4860  ) 

The  negative-ideal  alternative  is: 

R~  =  (  .3226.   .3145.   .3268.   .3913.   .4710.   .1608.   .1747,   .0694  ) 

££e.E   4^      Calculation    of    Separation    Measures 

Forthe  p=l  case, expressions  (2.5)  and  (2.7)  are  used  to 
calculate  d  (i)  and  d  "(i)  respectively,  and  the  result  is  given 
in   Tab  1  e    2.9. 

Table    2.9    Separation   Measures    When    p    =    1 
i  dj*(i)  dj"(i) 

1  .4584  1.2392 

2  .7365  .9611 

3  .5078  1.1888 

4  1.0142  .6834 

5  1  .  1866  .5110 

Similarly,  separation  measures  for  the  p  =  2  case  are  obtained 
and  are  shown  in  Table  2.10. 
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Table  2.10  Separation  Measures  When  p  =  2 
i       d2*U)  d2~(i) 


1 

.  2125 

2 

.4239 

3 

.  2455 

4 

.  4797 

5 

.5368 

.5961 
.4855 
.5833 
.3560 
.2659 


For  the  p  =  oo  case,  expressions  (2.6)  and  (2.8)  are  used  to 
calculate  d  *(i)  and  d  "  respectively,  and  the  result  is  shown  in 
Table  2.11. 


Table  2.11  Separation  Measures  When  p  =  oo 
i         d   *(i>  d   "(I) 


.  1165 
.  3496 
.1388 
.3217 
.4166 


.4166 
.4166 
.3496 
.  2778 
.2145 


Step_  5_;_  Unification  of  Two  Separation  Measures 

Unification  of  d_    is  obtained  using  (2.11).   For  eight 


attributes, 


.6154, 


'1  "  .«»■»»••   ^2 
Similarly.   unification  of  d 
results  are  given  in  Table  2.12. 


.2479.   and  C, 


.1367  are  used. 


is  obtained  using  (2.11).   The 
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Table  2.12  Unified  d    &  d    and  Their  Rankings 


d  ( i  )    ranks 


d  ( i )   ranks 


3536 

1 

6115 

3 

3956 

2 

7938 

4 

9282 

5 

9757 

1 

7756 

3 

9318 

2 

5516 

4 

4133 

5 

§.te£>  6j_  Ran.kj.ng  The  Preference  Order 

The  values  for  all  D's  and  the  final  ranking  order  are  shown  in 
Table  2.13. 

Table  2.13  Final  Solution 
i       D(i)  Rank 


*3' 

(Notice  that  Step  6  here  is  redundant,  because  a  unique  ranking 

order  has  been  achieved  at  Step  5.) 


1 

.  7340 

2 

.  5592 

3 

.7020 

4 

.4100 

5 

.3081 

The  sol 

uti 

on  is 

A, 

/ 


36 


2.3.3  Selection  of  Automobile 

A  customer  wants  to  purchase  a  family  car  and  is  faced  with 
the  problem  of  choosing  among  eight  automobiles.  The  ratings  based 
on  sixteen  attributes  are  given  in  Table  2.14.  The  data  is  adopted 
from  Fleischer  [1984]. 

The  given  problem  is  solved  by  the  following  steps: 

i.i£E  lx   Normalization  of  Decision  Matrix 

Expression  (2.5)  is  used  for  the  q  ■  1,  and  2  cases.  The 
normalized  matrices  are  shown  in  Tables  2.15  and  2.16.  Expression 
(2.6)  is  used  for  the  q  •  oo  case.  The  normalized  matrix  is  given 
in  Table  2.17. 

Stej>  2^   USifiSStiP.!!  °L   Xllll££  Normalized  Matrices 

The  coefficients,  Cp.  (  relative  reliabilities  )  for 
combining  matrices  are  determined  by  the  number  of  alternatives. 
Hence,  the  unified  normalized  matrix   is  obtained  using: 

Rli  -  .6154  rIj(l)  +  .2479  rij(2)  +  .1367  r^loo). 

Step.  3^   Unified  Decision  Matrix 

The  matrix  calculated  in  Step  2  is  shown  in  Table  2.18. 

§.l£Eij.  B£i£Efii!15lioil  °1  Positive- .Ideal,  and  Negative  -ideai 
Soiutions 
R*  =  (  .1231.  .1928.  .4442,  .3455,  .3798.  .3684.  .1988.  .2394. 
.3629,   .3391.   .3486,   .3675,   .3307.   .3247.   .3238.   .3546  I 

R"  ■  (  .4421.   .3295.   .1110.  .1234,   .1799,   .1973,   .3515.   .3222, 
.1561.  .0484.  .1307.  .1135.  .2480.  .2680,  .1835,  .1773  ) 
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Table  2.13  !ton»ali:ed  flatris  For  q  =  1 

0.068?  0.OB63  0.0633  0.1149  0.1782  0.1656  0.1264  0.1329  0.0752  0.1591  0.1429  0.0554  0.1111  0.1179  0.1071  0.0833 
0.0802  0.1136  0.0759  0.1149  0.1188  0.1210  0.1642  0.1218  0.1503  0.1591  0.1429  0.1142  0.1111  0.1429  0.0107  0.1250 
0.0877  0.1320  0.0759  0,1609  0.0990  0.1019  0.0929  0.1149  0.1503  0.1591  0.1633  0.1257  0.1481  0.1246  0.1429  0.1250 
0.1136  0.1460  0.1266  0.0575  0.0891  0.1146  0.1118  0.1045  0.0997  0.1591  0.1224  0.1793  0.1111  0.1179  0.1071  0.1250 
0.2475  0.1475  0.2532  0.1379  0.1089  0.1019  0.1178  0.1314  0.1748  0.0909  0.1020  0.1195  0.1481  0.1179  0.1429  0.1250 
M915  0.1136  0.075?  0.1609  0.1881  0.1733  0.1393  0.1255  0.0997  0.0909  0.1224  0.0953  0.1111  0.1429  0.1429  0.1667 
0.0B13  0.1136  0.0759  0.1149  0.1188  0.1210  0.1316  0.1218  0.1503  0.1591  0.1429  0.1406  0.1111  0.1179  0.107!  0.1250 
0.2292  0.1475  0.2532  0.1379  0.0990  0.0955  0.1161  0.1407  0.0997  0.0227  0.0612  0.1696  0.1481  0.1179  0.1429  0.1250 


Table  2.16  Normalized  Matrix  For  q  =  2 

0.171?  0.2408  0.152?  0.3165  0.4865  0.4569  0.3532  0.3769  0.2057  0.4217  0.3933  0.1501  0.3112  0.3333  0.3000  0.2326 
0.2001  0.3169  0.1835  0.3165  0.3243  0.3339  0.4589  0.3455  0.4115  0.4217  0.3933  0.3097  0.3112  0.4038  0.3000  0.3488 
0.2189  0.3684  0.1835  0.4430  0.2703  0.2812  0.2595  0.3259  0.4115  0.4217  0.4494  0.3407  0.414?  0.3521  0.4000  0.3488 
0.2834  0.4075  0.3059  0.1582  0.2432  0.3163  0.3123  0.2965  0.2727  0.4217  0.3371  0.4861  0.3112  0.3333  0.3000  0.3483 
0.6177  0.4116  0.6117  0.3797  0.2973  0.2812  0.3292  0.3728  0.4735  0.2410  0.2309  0.3240  0.4149  0.3333  0.4000  0,3488 
0.2284  0.3169  0.1835  0.4430  0.5135  0.4921  0.3892  0.3560  0.2727  0.2410  0.3371  0.2597  0.3112  0.4033  0.4000  0.4651 
0.2030  0.3169  0.1835  0.3165  0.3243  0.3339  0.3676  0.3455  0.4115  0.4217  0.3933  0.3812  0.3112  0.3333  0.3000  0.3488 
0.5721  0.4116  0.6117  0.3797  0.2703  0.2636  0.3244  0.3991  0.2727  0.0602  0.1685  0.4599  0.4149  0.3333  0.4000  0.3488 


Table  2.17  Noruaiiied  Matrix  For  q  ■  oo 

0.27B3  0.5850  0.2500  0.7143  0.9474  0.9286  0.7696  0.9444  0.4300  1.0000  0.8750  0.3088  0.7500  0.3256  0.750O  0.5000 
0.323?  0.7700  0.3000  0.7143  0.6316  0.6786  1.0000  0.B657  0.8600  1.0000  0.8750  0.6373  0.7500  1.0000  0.7500  0.7500 
0.3544  0.3950  0.3000  1.0000  0.5263  0.5714  0.5654  0.3164  0.8600  1.0000  1.0000  0.7010  1.0000  0.8721  1.0000  0.7500 
0.4589  0.9900  0.5000  0.3571  0.4737  0.642?  0.6806  0.7429  0.5700  1.0000  0.7500  1,0000  0.7500  0.8256  0.7500  0.7500 
1.0000  1.0000  1.0000  0.8571  0.5739  0.5714  0.7173  0.9339  1.0000  0.5714  0.6250  0.6667  1.0000  0.8256  1.0000  0.7500 
0.369B  0.7700  0.3000  1.0000  1.0000  1.0000  0.3482  0.8919  0.5700  0.5714  0.7500  0.5343  0.7500  1.0000  1.0000  1.0000 
0.32B6  0.7700  0.3000  0.7143  0.6316  0.6736  0.3010  0.8657  0.8600  1.0000  0.8750  0.7343  0.7500  0.3256  0.7500  0.7500 
0.9262  1.0000  1.0000  0.8571  0.5263  0.5357  0.7063  1.0000  0.5700  0.1429  0.3750  0.9461  1.0000  0.8256  1.0000  0.7500 


Table  2.18  Unified  Norssliied  Matrix 

0.1231  0.1928  0.1110  0.2463  0.3598  0.3421  0.2705  0.3043  0.1561  0.3391  0.3051  0.1135  0.2480  0.2680  0.2428  0.1773 
0.1432  0.2537  0.1332  0.2463  0.2398  0.2500  0.3515  0.278?  0.3121  0.3391  0.3051  0.2342  0.2480  0.3247  0.1835  0.2659 
0.1567  0.2949  0.1332  0.3455  0.1999  0.2105  0.1988  0.2631  0.3121  0.3391  0.3486  0.2576  0.3307  0.2832  0.3238  0.2659 
0.2029  0.3262  0.2221  0.1234  0.1799  0.2368  0.2393  0.2394  0.2069  0.3391  0.2614  0.3675  0.2430  0.26BO  0.2428  0,2659 
0.4421  0.3295  0.4442  0.2962  0.2199  0.2105  0.2522  0.3009  0.3629  0.193B  0.2173  0.2450  0.3307  0.2680  0.323B  0.2659 
0.1635  0.2537  0.1332  0.3455  0.3798  0.3684  0.29B2  0.2374  0.2069  0.193B  0.2614  0.1964  0.2480  0.3247  0.3238  0.3546 
0.1453  0,2537  0.1332  0.2468  0.2398  0.2500  0.2816  0.2789  0.3121  0.3391  0.3051  0.2882  0.2430  0.2680  0.242B  0.2659 
0.4095  0.3295  0.4442  0.2962  0.1999  0.1973  0.2485  0.3222  0.2069  0.0484  0.1307  0,3477  0.3307  0.2680  0.323B  0.2659 
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Stej)  5^   Calculation  of  Deviations  From  ideals 

When  distance  function  parameter  p  is  1.  the  distances  to  the 
pos i t i ve- i dea  1  and  the  negat i ve- i dea  1  are  summarized  in  Table 
2.19. 

Table  2.19  Separation  Measurements  for  p  =  1 
i  d1*(i)  dj'li) 

1  1.2435  1.5273 

2  0.9340  1.5633 

3  0.7802  1.9449 

4  1 . 0742  1 . 4625 

5  0.3405  1.5622 

6  0.7046  1.8420 

7  0.9452  1.6878 

8  0.6745  1.2583 

Similarly,  we  obtain  distance  measurements  for  the  p  «  2  and  oo 
cases.  The  results  are  displayed  in  Tables  2.20  and  2.21 
respectively . 

Step_   6_^    Combination    of    Three    Distance    Measurements 

The  coefficients  (relative  reliabilities)  are  determined  by 
the  number  of  attributes.  Therefore,  the  unified  distance 
measurements    are    obtained    using: 

d    =     .7121    dj     +     .2042    d9    +     .0837    d0o 

The    results    are    summarized    in    Table    2.22. 
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Table  2.20  Separation  Measurements  for  p  =  2 


d,  (i) 


d2-(i) 


0 .5370 
0 .4750 
0.4363 
0.  4879 
0.4942 
0.4510 
0.4335 
0.5847 


0 .5607 
0 .5328 
0. 6096 
0 . 5203 
0.5254 
0. 5554 
0.5513 
0.4947 


Table  2.21  Separation  Measurements  for  p  =  oo 


0.3332 
0.3110 
0.3110 
0.2221 
0 . 3190 
0.3110 
0.3110 
0 . 2907 


0.3190 
0 . 2989 
0 . 2907 
0. 2907 
0.3332 
0. 2786 
0 . 2968 
0.3332 
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Table  2.22  Unification  of  Distance  Measurements 
1  d  ( i )    ranks      d~ ( i )      ranks 

1 
2 
3 
4 
5 
6 
7 
8 


1 .0230 

8 

1 . 2288 

6 

0.7881 

6 

1 .2470 

5 

0.6707 

4 

1 .5338 

1 

0.8832 

7 

1 . 1720 

7 

0 . 3701 

1 

1 .2476 

4 

0 .6199 

2 

1 .4484 

2 

0 . 7876 

5 

1 . 3393 

3 

0  .  6240 

3 

1  .  0249 

8 

Step_  1-      EsaiSiSS  l!ie  Preference  Order 

Notice  that  the  ranking  order  for  d  is  different  from  the  one 
for  d~.  This  conflict  is  resolved  using  expression  (2.9).  The 
values  for  all  D's  and  the  final  ranking  order  is  given  in  Table 
2.23. 

Table  2.23  Final  Solution 
i  D ( i )    ranks 


0 

.5460 

8 

0 

.6127 

6 

0 

.  6958 

3 

0 

.  5703 

7 

0 

.7712 

1 

0 

.7003 

2 

0 

.6297 

4 

0 

.6215 

5 
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It  is  clear  from  Table  2.23  that  the  preferred  ranking  order  is 
Ac.  AR,  Ao,  A7,  Afl .  A5 
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2.4   Remarks 

(1)  TOPSIS  assumes  that  each  attribute  takes  either 
monotonica 1 1 y  increasing  or  decreasing  utility.  This  monotonici ty 
requirement  is  a  very  reasonable  assumption  [Keeney  and  Raiffa. 
1976].  A  nonmonotonic  utility  (e.g..  square  utility)  is  a  rare 
case,  such  as  the  optimum  number  of  rooms  in  a  house  or  the  blood 
sugar  count  in  a  human  body,  where  the  best  utility  is  located 
somewhere  in  the  middle  of  an  attribute  range. 

(2)  TOPSIS  considers  the  distance  from  both  the  positive- 
ideal  and  the  negative-ideal  points  simultaneously  by  taking  the 
relative  closeness  to  the  positive-ideal.  Thus,  the  DM  is  freed 
from  situations  where  whether  taking  the  positive-ideal  or  the 
negative-ideal  as  reference  point  is  uncertain.  This  feature  is 
particularly  appealing  to  the  DM. 

(3)  The  family  of  TOPSIS  solutions  may  be  difficult  to 
illustrate  graphically,  but  the  behavior  of  choosing  among  the 
solutions  can  be  explained  by  way  of  Indifference  (value)  curves 
[Keeney  and  Raiffa.  1976].  A  DM  is  assumed  to  give  equal 
preference  to  any  alternatives  located  on  the  same  indifference 
curve . 

The  typical  TOPSIS  indifference  curves  for  d2  and  d0o  are 
shown  in  Fig.  2.3  and  2.4  respectively.  Any  curves  with  Cj»  >  0.5 
are  convex  to  the  preference  origin,  which  indicates  the  property 
of  the  diminishing  MRS  observed  in  most  indifference  curves.  On 
the  other  hand,  indifference  curves  with  Cj*  <  0.5  are  concave  to 
the  preference  origin.  This  is  an  unusual  case  of  indifference 
curves,   but  it  may  be  interpreted  as  a  risk-prone  attitude 
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c=.? 


C=.«     c=.5 


Figure  2.3  Indifference  Curves  in  TOcSIS  with  d. 
Distance . 
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c=.9 


Figure  2.4  Indifference  Curves  in  TOPSIS  with  d 
Distance 
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resulting  from  a  pessimistic  situation. 

(4)  When  we  normalize  attribute  ratings,  the  parameter  q 
affects  the  degree  of  mutual  interdependency  among  alternatives 
for  given  set  of  attributes.  This  is  especially  true  when  q  «  oo. 
a  normalized  attribute  rating  is  independent  of  the  others  except 
for  the  best  rating.  The  lack  of  considering  other  alternative 
ratings  leads  to  the  lowest  reliability  of  d0Q. 

When  we  measure  the  distance  of  two  alternatives,  parameter  p 
plays  the  role  of  mutual  interdependency  among  attributes  which 
form  an  alternative.  When  p  ■  oo.  the  distance  between  two 
alternatives  is  independent  of  the  other  attributes  except  the  one 
where  the  two  alternatives  have  the  largest  deviation  from  one 
another.   This  leads  to  the  lowest  reliability  of  dQ0.  too. 

(5)  TOPSIS  solutions  depend  on  the  weight  vector  assigned  by 
the  DM  or  calculated  by  the  analyst.  Final  solutions  can  be  very 
different.  For  example.  the  final  solution  of  the  fighter 
aircraft  selection  problem  with  w=  (0.30.  0.20,  0.10.  0.10.  0.20) 
is  (A3.  A2.  Al.  A4).  while  the  solution  for  w=(1.0,  1.0.  1.0,  1.0. 
1.0,  1.0)  is  (A3.  Al.  A4,  A  2 )  . 

(6)  Although  the  TOPSIS  algorithm  requires  tedious 
calculations,  the  understandable  logic,  ease  of  programming,  and 
ease  of  computerization  make  it  one  of  the  best  MADM  methods. 
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Chapter    3         TOPSIS    -    MODM 
3. 1  Introduction 

Hwang  and  Yoon  [1981]  developed  the  Technique  for  Order 
Preference  by  Similarity  to  Ideal  Solution  (TOPSIS)  for  Solving 
Multiple  Attribute  Decision  Making  (MADM)  problems.  The  concept 
of  TOPSIS  was  then  adopted  to  formulate  a  new  method  for  solving 
Multiple  Objective  Decision  Making  (MODM)  problems.  This  new 
method  called  TOPSIS-MODM  tries  to  find  a  solution  that  is  closest 
to  the  positive  ideal  solution  (PIS)  and  farthest  to  the  negative 
ideal  solution  (NIS)  by  minimizing  some  relative  distance 
measurements . 

ESi^ti ve    Measure    of   Distance: 

To  measure  the  composite  distance  between  two  points  in  n- 
dimensions  is  almost  impossible.  because  we  are  facing  the  issue 
of  noncommensurabi 1 1 ty  of  Individual  criteria.  It  is  necessary  to 
make  all  dimensions  commensurable  and  calculate  the  relative 
deviation    instead    of    absolute    deviation    between    two    points.        The 

1  2 

function  of  the  relative  measure  of  distance  between  f_   and  f_      can 
be  written  as : 


-  f. 


%  -  {  =  wr 


f      }1/P 


(3.1) 


where   Wr    is    the    relative   weight    assigned    by    the   DM,    and   p   can   be 
any    integer    between    1    and    oo    . 
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For  p  =  oo  ,  the  above  expression  is  rewritten  as: 


.    1    _  t    2 
rr     lr 


fr   -  f, 


(3.2) 


The  property  of  p  is  best  described  by  Yu  [1985]:  "...  p  plays 
the  role  of  the  balancing  factor  between  the  "group  utility  sum" 
and  the  maximum  of  the  individual  regrets.  If  we  want  to 
emphasize  the  utility  sum  we  should  select  a  small  p  ;  otherwise, 
we  should  select  large  p  in  order  to  increase  the  weight  to  the 
individual  regrets...."  Boychuk  and  Ovchinnikov  [1973]  have 
suggested  p  ■  1,  and  Salukvadze  [1972]  has  suggested  p  =  2.  The 
most  widely  used  p.  however,  are  1.2.  and  oo.  (See  Yoon  [1986], 
Zeleny  [1982]  and  Cohon  [1978]). 

The  function  of  relative  measures  of  distance  are  used  later 
in  this  chapter  to  measure  two  distances:  (1)  the  distance  to 
the  positive  ideal  solution,  and  (2)  the  distance  to  the  negative 
ideal  solution. 

TOPSIS-MODM  is  similar  to  the  global  criteria  method  in  the 
sense  of  finding  the  solution  that  is  closest  to  the  positive 
ideal  solution.  However,  the  global  criteria  method  did  not 
consider  that  a  solution  which  is  farthest  to  the  negative  ideal 
solution  can  also  be  a  preferred  solution.  The  composite  of  these 
two  distance  measurements,  i.e.  distance  to  the  positive  ideal  and 
the  distance  to  the  negative  ideal,  are  the  main  features  of  this 
new  method. 
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Formulation   of   TOPSIS-MODMj. 


(i)    Determine    the    alternative    that    is    closest    to    the    PIS: 


fj       "    fj(S)       p  'i<£)     "    *i  p 

lin    F_    =    {    Z    w.     [ . ]       +       IHj     [ 1       > 

V  )        J  -  , 


s.  t.     x  €  X  (3.3) 

(ii)  Determine  the  alternative  that  is  farthest  to  the  NIS: 

max  i,  .  {  i  «,  [    J        __i r    +    s  w,  [„ r  > 


fj    -  fj"  fi"  ~  fi 


,t.     x    £    X.  (3.4) 


The       combination    of    equations    (3.3)    and    (3.4)     for    p    -    1.2 oo 

cases    can    be    generalized    as: 

Kin    (Hp    Fp    -    Hn    Np)  (3.5) 

s. t .  x    €    X 

Where  H  and  Hn  are  the  importance  indexes  of  distances  to  PIS  and 
NIS.  respecti ve  1  y.  (see  next  sectilon).  However,  when  formulating 
(3.5)  for  the  p  »  oo  case,  we  find  that  it  is  not  possible  to 
combine  (3.3)  and  (3.4)  without  transforming  (3.3)  to  (3.6)  and 
(3.4)  to  (3.7)  respectively. 
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min  d 


(3.6) 


,t.   x  6  X 


WJ  t 


"it. 


and 


fj    -  fj(x) 


fj   "  f i 


fi<x)  -  fj 

* 

f.    -  fi 


<  d      (All  j  6  J) 


.]  <d      (All  16  I! 


max  d" 


(3.7) 


s.t.    x  €  X 


fj(*>  "  fj 


fj   -  fJ 


]  >  d" 


(All   j  6  J ) 


■  i   [ 


]  >  d" 


(All   i  €  I) 


And  merging  (3.6)  with  (3.7).  we  obtain  the  formulation  for  p 
oo.  (see  expression  (3.15)  for  details.) 
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Introducing  the  Importance  Index  H_  and   Hn: 

As  mentioned,  TOPSIS-MODM  tries  to  find  a  solution  that  is 
closest  to  the  PIS  and  farthest  from  the  NIS.  Unfortunately,  such 
a  solution  may  never  exist.  In  order  to  obtain  a  satisfactory 
solution,  a  set  of  importance  index  H  (H  ,  Hn)  are  incorporated  in 
to  the  TOPSIS-MODM  formulation.  Hp  and  Hn  have  the  following 
properties : 

(1)  H  (Hn)  is  an  importance  index  assigned  by  the  DM  to 
indicate  the  importance  of  choosing  a  solution  that  is  closest 
(farthest)  to  the  PIS  (NIS). 

(2)  H  and  H„  are  normalized  weights,  hence,  their  values 
will  always  be  between  0  and  1,  i.e.  0  <  Hp  <  1 .  0  <  Hn  <  1  .  and 
the  summation  of  these  two  indexes  should  always  equal  to  1,  i.e. 

"p  +  Hn  "  »■ 

(3)  The  values  of  H  and  Hn  may  not  be  the  same.  Larger  H„ 
implies  the  DM  wants  to  use  the  PIS  as  the  reference  point  while 
smaller  H„  implies  the  DM  favors  the  NIS  to  be  the  reference 
point.  If  H  equals  H  ,  then  we  conclude  that  a  compromise 
solution  referring  to  both  the  PIS  and  the  NIS  is  desired  by  DM. 

(4)  When  H_  •  1  and  Hn  =  0,  TOPSIS-MODM  is  the  method  of 
global  criteria.  In  other  words,  the  method  of  global  criteria  is 
a  special  case  of  TOPSIS-MODM. 

(5)  Although  H  and  Hn  values  have  significant  effects  on 
TOPSIS-MODM,  one  exception  is  observed.  When  p  =  1  is  used  in  the 
TOPSIS-MODM  formulation,  H  and  Hn  values  have  no  effect  on  the 
final  solution,  i.e.  H  =  (0.75.  0.25)  and  H  -  (0.50.  0.50)  give 
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identical  solutions  for  the  p  =  1   case  (See  the  first  numerical 
example  for  details.) 
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3.2   Algorithm  for  TOPSIS-MODM 
P£.°kle™  E°£!!l*il  Stion^ 

Let  the  MODM  problem  be  : 

Max/  Min  [  f  j (x) .  f 2(x) . 

s.t.   x  €  X 


fk(x)] 


(3.8) 


where  X  =  {  x  |  gj(x)  {>.=.<}  0.  i  =  1.  2 m  } 

fj(x):  the  benefit  objective  for  maximization,  j  € 
fj(x):  the  cost  objective  for  minimization.  1  €  I 


The  complete  TOPSIS-MODM  algorithm  for  this  problem  is  given 
below . 

Step.  <K   Determine  the  Positive  Ideal  Solution  (PIS),  f  .  and 
the  Negative  Ideal  Solution  (NIS).  f". 

Step.  0-l_^  Obtain  and  Construct  PIS  -  Payoff  Table 
The  problem  is  to  find  x  so  as  to: 


max  fj(x)  .  for  all  j 


s.t.    x  6  X 
and  /  or 

min  f  j  (x) ,  for  all  i 
s.t.   x  €  X 


(3.9) 


(3.10) 


The  solutions  of  (3.9)  and  (3.10).  f j   and  fj   are  known  as 
the  positive  ideal  solutions.   Each  positive  ideal  solution,  e.g. 
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this  set  of  optimal  decision  variables, 
objective  functions  fj,   f2 'r-1"   ^r  +  1 


into  the  other 


we  obtain 


'r.l-   fr.2 fr.r-l'   fr,r  +  l fr.k-   Then  the  PIS  P"?"" 

table  is  constructed  as  follows: 

PIS  -  Payoff  Table 


l21 


12 


Ik 


2k 


•kl 


_k2 
to' 


(  f 


Steg    0-2^      Obtain    and    Construct    NIS-Payoff    Table 
The   problem   is   to   find   x   so   as   to 


min    f  -(x) .     for    all     j 
s. t.         x    6    X 


(3.11) 


and    /    or 


max    f i (x) ,      for    all     i 


(3.12) 


s.t. 


x    €    X 


The    solution    of    (3.11)    and    (3.12).    x"   and   f, 
negative    ideal    solutions. 


ire    known    as    the 


( 
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With   all    the   x      from    (3.11)    and    (3.12).    the   NIS-Payoff   table    is 
gi  ven    as : 

NIS-Payoff   Table 


'21 


'12 


Ik 


'2k 


'kl  zk2 


r  -  (  ff.  f2  'k    ' 

*         _ 

Note  that  f  and  f  provide  the  upper-  and  lower-bounds  on  the 
objective  function  values  and  are  then  important  reference  values 
for  the  DM  in  all  the  following  steps. 

Step  1.    Calculation  phase: 

The  optimal  solution,  x  and  f.  can  be  obtained  by  solving  the 
following  problems. 

(1)   When  p  »  1 .   The  problem  is  to  find  x  so  as  to: 
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min    {    H       [    S    Wi 
i 


fj       -    fj(x) 


fj       "    fJ 


fjlx)    -    f 


*ll 


Hn    t?    w.i 


fj<*>    "    fj" 


fj 


fi"   "    'l<D 


]    } 


.  t.         x    6    X 


(3.13) 


(2)       When    p    ■    2.       The    problem    is    to    find    x    so    as    to: 


»in    {    H       [    S   W1 


fj       -    fj(x) 


fj       -    f i 


Z    W( 


_  * 

t,       -    t. 


Hn    I    ?    wj 


fj<*>    "    h 


fj" 


+       E    W, 


fi        "    fi 


]     > 


.t.  x    €    X 


(3.14) 


(3)      When   p   =   oo    .      The   problem    is    to    find   x    so    as    to: 

min     (     d*     -    d~)  (3.15) 

s.t.  x    €    X 
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where 


Hp  Wj  [___ ]     <  d    (  j  €  J  ) 

fi(*>  "  fi* 
Hp  W j  [ ]     <  d   (  i  6  I  ) 

Hn  Wj  [__• i ]    >  d"  (  j  e  j  ) 

fi"  -  fi<x) 
Hn  Wj  [ ]     >  d   (  i  €  I  ) 

_      * 

(a)  Wj.   i  «  1.2 k  are  weights  (importances)  of 

the  objectives. 

(b)  H   is  the  DM's  weight  on  how  important  it  is  to 

obtain  a  solution  that  is  closest  to  the  PIS.  Hn  is 
the  DM's  weight  on  how  important  it  is   to  obtain  a 
solution  that  is  farthest  from  the  NIS. 
(C)    Hp  +  H„  =  1 


Step  2.    Decision  phase. 

The  solutions  obtained  from  previous  step  are  presented  to  the 
DM.   The  DM  may  select  one  or  none  of  them. 
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3.3  Numerical  examples 

This  section  contains  three  numerical  examples.  The  first  one 
is  a  two- var i ab 1 e  .  two -ob j ec t i ve  linear  model  which  is  used  to 
graphically  demonstrate  the  power  of  TOPSIS-MODM.  The  second  one 
is  a  much  more  complex  linear  problem  which  is  more  realistic  than 
the  first  example.  The  third  one  is  a  nonlinear  model  to  show  that 
TOPSIS-MODM  is  also  good  for  nonlinear  cases. 

3.3.1.  Graphic  Example:  Hardee  Company  Production  Scheduling 

The  Hardee  company  makes  two  kinds  of  dolls.  Doll  A  is  a  high 
quality  one  and  Doll  B  is  of  lower  quality.  The  respective 
profits  are  $0.40  and  $0.30  per  doll.  Each  Doll  A  requires  twice 
as  much  time  as  a  Doll  B.  and  if  all  dolls  were  of  type  B.  the 
company  could  make  500  per  day.  The  supply  of  material  is 
sufficient  for  only  400  dolls  per  day  (both  A  and  B  combined). 
The  problem  assumes  that  all  the  dolls  for  type  A  and  type  B  the 
factory  can  make  could  be  sold,  and  that  the  best  customer  of  the 
company  wishes  to  have  as  many  as  possible  of  type  A  doll.  The 
manager  realizes  that  two  objectives:  (1)  the  maximization  of 
profit,  and  (2)  the  maximum  production  of  Doll  A.  should  be 
considered  in  scheduling  the  production. 

This  M0DM  problem  can  be  formulated  as  follows: 


max  f2(x)  =  Xj^ 
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subject  to 


g,  ( x)  =  x,  +  x2  <  400 


g,(x)  =  2  x,  +  x,  <  500 


x j .  X2  >  0 


produced  respectively.   And  let  constraints  be  denoted  by   x  €  X. 

The  steps  of  the  TOPSIS-MODM  algorithm  are  as  follows. 

Si£E  2j.  2klSi!l  afi<i  £°fistruct  Pay_of  f  Tables 

0.1  (PIS)     1)  max  f](x)  =  0.4x:  +  0.3x2 

s.t.   x  €  X 

The  solution  is:  Xj  =  100.  x2  =  300;  fj(x  )  =  130 
Point  B  is  the  ideal  solution  for  f j  (x)  (See  Figs.  3.1  and 
3.2)  . 

Similarly,  for  the  second  objective: 

2  )  max  f  j ( x )  «  Xj 

s.t.   x  €  X 

The  solution  is:  Xj  »  250.  x2  ■  0:  f2(x  )  =  250 
Point  C  is  the  ideal  solution  for  f2(x)  (See  Figs.  3.1  and  3.2). 
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The  PIS-payoff  table  is  : 


1 


f; 


130 
100 


100 
250 


100 
250 


300 

0 


Note  that  in  this  payoff  table,  row  1  gives  the  solution 
vector  x  =  (100.  300)  which  maximizes  the  objective  fi(x  )  =  130; 
f2  (100,  300)  =  100  is  the  value  taken  on  by  the  objective  f2  when 
f j  reaches  its  maximum.  Similarly,  row  2  gives  the  solution 
vector  x  »  (250,  0)  which  maximizes  f2(x  )  ■  250;  f,(250,  0)  • 
100  is  the  value  taken  on  by  f j  when  f2  reaches  its  maximum.   The 


point  P .   (  f  j 


)  =  (130.   250).   in  Figure  3.2  is  an  Utopia 


solution  called  the  PIS.  which  is  not  feasible. 
0-2.  (NIS)       1)  min  fjlx)  =  0.4x1  +  0.3x2 

s.t.  x  e  x 


The  solution  is: 


■1 


s.t.  x  e  x 
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The  NIS-payoff  table  is: 


fl 

f2 

xl 

x2 

'l 

0 

0 

0 

0 

f2 

120 

0 

0 

0 

<    *2 

< 

400 

The  NIS  is  point  N,  (fj  .  f2  )  ■  (0.  0).  in  Figure  3.2. 


Step,  Ll  SkiSifi  the  p_referred  solution 

A  preferred  solution  is  a  nondomina t ed  solution  which  is  a 
point  on  the  segment  of  straight  line  BC.  (In  Figs.  3.1  and  3.2, 
all  points  on  the  segment  of  straight  line  BC  are  nondom i nated 
solutions).   Assume  equal  weights  for  the  two  objectives  and  H   ■ 


0.75  and  H 


n 


0.25. 


Case  1 .  P  =  1 

The  problem  becomes  to  find  x  so  as  to: 


130  -  0.4x, 


0.3xo 


min  {  (  0.75  )  [- 


130 


250  -  x, 


(  0.75  )  [ ] 

250  -  0 


(  0.25  )   [• 


130  -  0 


(  0.25  )   [ ]     ) 

250  -  0 
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Figure  3.1  TOPSIS-MODM  for  Hardee  Co.  Production  Scheduling 
Problem  with  H  =  (0.75,  0.25),  W  -  (0.50,  0.50) 


t1or&o«\inanA     £,/*  lions 


ICO    no 


Figure  3.2  Objective  Function  Representation  of  Feasible 

Area  of  Hardee  Co.,  Production  Scheduling  Problem 


s.t.    x  €  X 

This  is  a  LP  problem.  The  solution  is  : 

x  =  (  250  ,  0  ) 

f  =  (  100,  250  ) 
which  is  point  C  in  Figs.  3.1  and  3.2. 

Case  2 .  P  =  2 

The  problem  becomes  to  find  x  so  as  to: 


nin  {  (  0.75  )  [■ 


+  (  0.75  )  [- 


130  -  0.4X,  -  0.3x 


2   ,2 


130  -  0 


250 


(  0.2  5  )   [■ 


(  0.25  )  [■ 


s.t.   x  €  X 


This  is  a  nonlinear  programming  problem.  We  use  Sequential 
Unconstrainted  Minimization  Technique  (SUMT)  to  solve  it.  Tlie 
solution  is 

X  ■  (  232.21.  35.61  ) 

f  =  (  103.56,  232.21  ) 
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which  is  point  D  in  Figs.  3.1  and  3.2. 
Case  3.  P  =    oo 

The  problem  is  to  find  x  so  as  to: 

*     _ 
min  {  d   -  d   } 

s.t.   x  6  X 

130  -  0.4x,  -  0.3x2 
(  0.75  )  [ ]    <  d 

130  -  0 

250  -  X, 
(  0.75  )   [ ]      <  d* 

250  -  0 

0.4x1  +  0.3x2  -  0 
(  0.25  )  [ ]      >  d" 

130  -  0 

Xj  -  0 
(  0.25  )  [ ]      >  d" 

250  -  0 

This  is  a  LP  problem.  The  solution  is  : 
X  =  (  208.74,  82.52  ) 
f  =  (  108.25,  208.74  ) 
which  is  point  E  in  Figs.  3.1  and  3.2. 

Stej)  2^  Decision  Phase 

The  solutions  for  p  =   1.2  and  oo  cases  are  summaried  in 
Tab le  3.1  .  ' 
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Table  3.1  Solutions  Table 


PIS 

130.00 

250. 00 

NIS 

0 

0 

1 

100 .00 

250.00 

250 

00 

0.00 

2 

104 .00 

230 . 00 

230 

00 

40.00 

00 

108.25 

208.74 

208 

74 

82.52 

The  solutions  table  is  presented  to  the  DM.  The  DM  may  select  one 
or  none  of  the  solutions. 


As  shown  in  Fig.  3.1  for  p  =  1,  the  function  to  be  minimized 
F„_j.  is  a  linear  function;  therefore  the  preferred  solution  is 
determined  by  the  slope  of  line  BC  and  that  of  F  j.  Since  the 
slope  of  F  j  is  less  than  that  of  line  BC  (the  line  representing 
the  nondominated  solutions),  point  C  is  the  preferred  solution. 
For  p  =  2,  the  function  to  be  minimized.  F  2'  *s  a  nonlinear 
(quadratic)  function,  and  the  minimum  point  for  F  2  Is  point  D  in 
Fig.  3.1.  which  is  the  point  where  the  F  2  curve  touches  the 
segment  of  line  BC.  For  p  =  oo.  the  function  to  be  minimized. 
F  00.  cannot  be  drawn  on  Fig.  3.1  but  can  be  represented  in  Fig. 
3.2  as  point  E. 


Note: 

Several  sets  of  compromise  solutions  are  obtained  in  Table 
3.1.1  for  the  Hardee  Company  production  problem  by  utilizing 
different  H   and  Hn  combinations. 
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Table  3.1.1   Solution  Tables  for  W  =  (0.50.  0.50) 


PIS 
NIS 


f  1 
130. 00 

0 

f  1 


f2 
250.00 

0 
f2 


xl 


x2 


xl 


x2 


(1  .00 ,  0.00) 


1 

2 
oo 


100.00  250.00  250.00  0.00 
103.90  230.70  230.70  38.60 
108.18     209.09      209.09     81.82 


(0.75,  0.25) 


1 

2 

00 


100.00  250.00  250.00  0.00 
103.56  232.21  232.21  35.61 
108.25     208.74       208.74     82.52 


(0.50.  0.50) 


100.00  250.00 
103.21  233.44 
107 .97     210.  14 


250.00       0.00 
233.44     32.78 

210.14     79 . 72 


(0.25 .  0.75) 


1 

100 

00 

250 

00 

250 

00 

0 

00 

2 

100 

.  21 

248 

.98 

248 

.98 

2 

.05 

00 

108 

.25 

208 

.  74 

208 

74 

82 

52 

1 

100 

00 

250 

.  00 

250 

00 

0 

00 

2 

100 

,00 

250 

.  00 

250 

,  00 

0 

.00 

00 

108 

18 

209 

09 

209 

0  9 

81  . 

32 

(0. 00.  1 .00) 


It    is    obvious    that    when    parameter    p    equals    1,     the    solutions 
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for  the  problem  are  always  the  same,  i.e.  H  has  no  effect  on  the 
problem's  solution  when  p  «  1,  That  Is.  the  slope  of  Fp»l  remains 
constant  and  is  less  than  that  of  line  BC;  therefore  the  preferred 
solutions  are  shown  at  the  same  point.  C.  in  Figs  3.1  and  3.2. 
Almost  all  the  solutions  for  p  =  2  and  oo  are  changed 
correspondingly  when  H  is  changed.  Notice  also  that  when  H  ■ 
(1.0,  0),  the  set  of  solutions  is  exactly  the  same  as  with  the 
method  of  global  criteria.  (See  the  first  numerical  example  of 
appendix    B.) 
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3.3.2.    Nutrition  Problem 

This  problem  is  to  determine  the  quantities  of  six  foods  that 
should  be  eaten  to  balance  nutritional  requirements  so  as  to: 
(i)  maximize  carbohydrate  intake 
(ii)  minimize  cholesterol  intake 
(iii)  minimize  cost 

Information  on  six  foods  is  given  in  Table  3.2. 

Table  3.2.   Nutrition  and  cost  information  on  six  foods 


Vitamin  A  ( i  .  u  .  ) 

Food  energy 
(calories  ) 


Lettuce   Orange   Daily 
Milk   Beef    Eggs    Bread   salad    juice   intake 

for 
(pint)  (lb.)  (dozen)  (oz.)   (oz.)   (pint)   adults 

720    107   7080      0      134      1000      5000 

344   1460   1040     75       17.4     240      2500 


Cholesterol 

10 

20 

120 

0 

0 

0 

(unit) 

Protein  (g) 

18 

151 

78 

2 

.5 

0 

.2 

4 

Carbohydrate ( g ) 

24 

27 

0 

15 

1 , 

1 

52 

Iron  (mg) 

0. 

2 

10 

1 

13 

2 

0 

75 

0 

15 

1  . 

2 

Cost  ($) 

0  . 

22 

2. 

2 

0. 

8 

o . 

1 

0 

05 

0. 

26 

12.5 
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The  problem  is  summarized  as: 

max  fj(x)  =  24Xj+  27x2  +  15x4  +  l.lx5+  52xg  (Carbohydrate  intake) 


min  f,(x)  ■  10x,+  20x9  +  120x, 


(Cholesterol  intake) 


min  f3(x)  =  0.22Xj+  2.2x2+  O.8X3+  "  ■  lx4+  °-05x5+  0.26x6  (Cost) 
Subject  to  : 

720xj  *  107x2  +  7O8OX3  +  134x5  +  lOOOXg  >  5000  (Vitamin  A.i.u.) 
0.2Xj  +  10.1x2  +  13.2x3  +  0.75x4  +  0.15x5+  1 . 2xg  >12.5  (Iron.mg) 
344Xj  +  1460x2  +  1040x3  +  75x4  +  17.4x5+  240x6  >  2500  (Calories) 


I8X1  +  151x,  +  78Xo  +  2.5X.,  +  0.2x=  +  4xK  >  63 


'1 

£1  1  6  •  ° 
:,  <  1.0 


6  i 


( Protein,  g) 


xfi  <  4.0 


(Milk,  pints) 
( Bef f  .  pounds  ) 
( Eggs .  dozen) 
(Bread,  ounces) 
(Lettuce  &    Sal  ad , ounces ) 
(Orange  juice,  pints) 


Xj  >  0  ,  i  =  1 6 

J  =  <1>.   I  ={2,  3) 
Let  the  constraints  set  be  denoted  by  x  6  X 


The  problem  is  solved  by  the  following  steps. 
Step  0.    Obtain  and  Construct  PIS-  &  N^S- Payoff  Tables 

LI.s-E2y.0ff  table 


1 


'1 


max  f  j  5  4  0 

00 

80 

00 

6 

09 

6 

00 

1  .00 

0  .00 

10 

00 

10 

00 

4  .00 

min  f2  380 

3  9 

8 

44 

3 

47 

0 

00 

0.42 

0.00 

10 

00 

10 

00 

4.00 

min  f3  281 

67 

67 

83 

2 

26 

3 

78 

0.00 

0  .  25 

10 

00 

0 

00 

0 .  79 
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f   =  (540.00.  8.44,  2.26) 
NJ.S-  P_ay_off  table 


max  f  j 

93 

34 

71 

53 

3.08 

2 

50 

0.83 

0 

25 

0 

00 

10 

00 

0 

00 

min  f o 

321 

00 

110 

00 

4.75 

6 

00 

1  .  00 

0 

25 

10 

00 

0 

00 

0 

00 

m  in  f  3 

540 

00 

110 

00 

6  .  29 

6 

00 

1.00 

0 

25 

10 

00 

10 

00 

4 

00 

f   =  (93.34 ,  110.00  .  6 . 29) 

ii£E  k—   Calculation  fihase 

We  assume  the  DM  equally  weights  the  objectives  and  assign 
Hp  -  0 .75  and  Hn  =  0. 25. 

Case  1 .  p  =  1 

The  problem  is  to  find  x  so  as  to: 

540.00  -  fjtx) 


«in  {  (  0.75  )  [- 


+  (  0.75  )  [- 


+  (  0.75  )  [■ 


540.00  -  93.34 


f,(x)  -  8.44 


110.00  -  8.44 


f, (x)  -  2 .26 


L3'£ 


6.29  -  2 . 26 
fj(x)  -  93.34 


-  (  0.25  )  [ ] 

540.00  -  93.34 

110.00  -  f,(x) 

-  (  0.25  )  [ ] 

110.00  -  8 .44 

6.29  -  f3(x) 

-  (  0.25  )  [ ]  > 

6.29  -  2.26 


s. t .   x  €  X 
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Since  all  the  objective  functions  and  constraints  are   linear,  the 
problem  is  a  linear  model.  The  solution  is 

X  -  (  2.30,  0.00,  0.00,  10.00,  0.00.  4.00  ) 

f  =  (  413. 12.  22.97.  2.55  )  . 

Case  2.  p  »  2 

The    problem    becomes    to    find    x    so    as    to: 

540.00    -    fjtx)     2 

min    {     (    0.75     )     [ ] 

540.00    -    93.34 

f 2(x)     -    8. 44  2 

+    (    0.75     )     [ ] 

110 .00    -    8 . 44 

f 3 (x)     -2.26  2 

*     (     0.75     )     [ ] 

6. 29    -    2 . 26 

f i (x)     -    93 .34       2 

-  (    0.25     )     [ ] 

540.00    -    93.34 

110.00    -    f2(x)    2 

-  (    0.25     )     [ ] 

110.00    -    8.44 

6.29    -    f,(x)  2 

-  (    0.25     )     [ ]"     } 

6. 29    -    2 . 26 

s  .  t .     x    S    X 

This   is   a  nonlinear  problem.   We   use  SUMT  to  solve   it.   The 
solution    is 

X    =     (     1.69.     0.09.     0.01,     9.99,     3.61.     4.00     ) 
f    -     (     404.83  .     19 . 68 .     2.81     )  . 
Case    3 .    p    =    oo 

The    problem    is    to    find    x    so    as    to: 
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min  (  d   -  d"  ) 

s. t .    x  €  X 

540.00  -  f. (x) 

(  0.75  )  [ ]      <  d 

540.00  -  93.34 

f2(x)  -  8.44 
(  0.75  )  [ ]       <  d* 

110 . 00  -  8. 44 

f3(x)  -  2  .  26 

(  0.75  )  [ ]       <  d 

6.29  -  2 .26 

fjlx)  -  93.34 
(  0.25  )  [ ]      >  d" 

540.00  -  93.34 

110. 00  -  f, (x) 

(  0.25  )  [ ]      >  d" 

110. 00  -  8.44 

6.29  -  f3(x) 
(  0.25  )  [ ]       >  d~ 

6. 29  -  2 . 26 

This  is  a  LP  problem.  The  solution  is 

X  =  (  3.08.  0.00.  0.00,  10.00.  8.68.  4.00  ) 
f  =  (  441  . 47.  30. 80  .  3.15  ) 

Steg    2^_   Decision    Phase 

The    solutions    for    p    ■    1,2    and    oo    cases       are    summarized    in 
Tab  1  e    3.3. 
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Table  3.3  Solutions  Table 


PIS 

540 

00 

8 

44 

2 

26 

--- 

--- 

... 

— 

— 

— 

NIS 

93 

34 

110 

00 

6 

29 

— 

... 

-— 

--- 

--- 

... 

P  -  1 

413 

12 

22 

97 

2 

55 

2  .30 

0.00 

0 

00 

10  .00 

0  .00 

4  .00 

2 

404 

83 

19 

68 

2 

81 

1  .69 

0  .  09 

0 

01 

9  .99 

3.61 

4  .00 

00 

441 

47 

30 

80 

3 

15 

3.  08 

0.00 

0 

00 

10.  00 

8  .  68 

4.00 

The  solutions  table  is  presented  to  the  DM  for  evaluation. 
The  DM  may  select  one  or  none  of  the  solutions. 
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3.3.3    Bow  River  Valley  Water  Quality  Management  Problem 

This  problem  is  based  on  a  hypothetical  case  developed  by 
Dorfman  and  Jacoby  [1969]  and  modified  by  Monarchi  et  al.  [1973]. 
It  is  concerned  with  the  pollution  problems  of  a  river  basin,  the 
Bow  River  Valley,  whose  main  features  are  shown  in  Fig.  3.3. 
Industrial  pollution  originates  from  the  Pierce-Hall  Cannery, 
located  near  the  head  of  the  valley,  and  municipal  waste  pollution 
comes  from  two  outlets  at  Bowville  and  Plympton.  A  state  park  is 
located  between  the  cities,  and  the  lower  end  of  the  valley  is  a 
part  of  the  state  boundary  line. 

The  Bow  Valley  Water  Pollution  Control  Commission  is  made  up 
of  representatives  from  all  three  waste  dischargers  and  members  of 
the  state  and  federal  government.  The  commission  is  responsible 
for  maintaining  river  water  quality  by  setting  the  waste  reduction 
requirements  at  the  three  sources  of  pollution,  but  it  must  act 
with  an  awareness  of  the  effect  of  any  additional  effluent 
treatment  costs  on  the  economy  of  the  valley. 

The  Commission  (the  DM)  is  concerned  with  the  following  six 
objectives:  the  DO  levels  at  Bowville.  Robin  State  Park,  and 
Plympton  (objectives  1,  2.  and  3.  respectively);  the  percent 
return  on  investment  at  the  Pierce-Hall  Cannery  (objective  4):  and 
the  addition  to  the  tax  rate  for  Bowville  and  Plympton  (objectives 
5  and  6,  respectively).  The  first  problem  of  the  DM  is  to 
determine  a  waste  treatment  policy  which  satisfies  the  constraint 
that  the  DO  level  at  the  state  line  is  above  3.5  mg/1. 

The  problem  is  summarized  as  follows:  determine  Xj,  x,  and  x3 
so  as  to  satisfy  the  following  objectives,  constraints  and  goals. 
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Bow  River 
0 


(4.75) 


Bowvllle 

Population 

250.000 


(5.1) 


Plympton 

Population 

200.000 


10 


50 


70 


100 


P1erce-Hall 
Cannery 


Robin 
State  Park 


STATE  LINE 


(6. 75) 


(2.0) 


(1.0) 


Distance  1n 
river  miles 


(Values   1n  parentheses  at  the  sides  are  current 
00  levels   1n  milligrams  per  liter.) 


Figure    3.3   Main   Features   of  the   Bow  River  Vallev 
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The  six  objectives  are: 

max  f j       «      4.75      +      2.27       (      Xj       -      0.3       )  (      mg/1      DO) 

max  f2    =    2.0    +    0.524     (    Xj    -    0.3)     +    2.79     (    x2    -    0.3) 

+    0.882     (    Wj    -    0.3    )     +    2.65     ( w2    -    0.3)     (    mg/1    DO) 

max  f3    =    5.1    +    0.177     (    Xj    -    0.3)     +    0.978     (    x2    -    0.3) 

+    0.216     (wj    -    0.3)     +    0.768     (    w2    -    0.3)     (    mg/1    DO) 

59 
max  f4    «    7.5    -    0.012    (    -    59       )     (    X    return) 

1  .  09  -  Xj2 

532 

min     f5  -  1.80  *  10"3  ( 532)(  *  additional 

tax  rate) 
1.09  -  x22 

450 

min     f6  ■  2.50  *  10"3  (  -  450)(  *  additional 

tax  rate) 
1  .  09  -  x32 

s.t.   gj  =  1.0  +  0.0332(xj  -  0.3)  +  0.0186(x2-  0.3) 
+  3.34(x3  -  0.3)  +  0.0204(Wj  -  0.3) 
+  0.778(w2  -  0.3)  +  2.62(w3  -  0.3)  >  3.5 

0.3  <  Xj  <  1.0.  i  =  1.2.3  (3.21) 

Let  the  above  constraints  be  denoted  by  x  6  X . 
w^  is  related  to  Xj  by: 

0  .39 
wt  = .  1  -  1,2.3 

1  . 39  -  Xj2 
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The  original  problem  formulation  is  slightly  modified  to 
reduce  its  complexity.  Both  objectives  corresponding  to  DO  levels 
at  Bowville  and  Plymton  are  changed  to  constraints  for  DO  > 
6.0.   The  modified  problem  is  given  below. 

max  fjlx)  -  2.0  +  0.524(Xj  -  0.3)  +  2.79(x2  -  0.3)      (3.22) 

+  0.882IWJ  -  0.3)  +  2.65(w2  -  0.3)(mg/l  DO.  at 

state  park) 

59 

max  f2(x)  =  7.5  -  0.012  ( -   59)(«  earning,  for 

cannery) 
1.09  -  Xj 

532 

mln  f3(x)  =  1.80  *  0.001  (  -  532)  (  $/$1000 

addi ti tional 
1.09  -  x22  tax  at 

Bowvil le ) 

450 

min  f4(x)  =  2.50  *  0.001  (  -  450)  ($/$1000 

addi  tional 
1.09  -  x32  tax  at 

Plymton) 

s.t.  gx(x)  =  4.75  +  2.27  (  Xj  -  0.3)  >  6.0 

g2(x)  =  5.1  +  0.177  (  Xj  -  0.3)  +  0.978(x2  -  0.3) 
+  0.216(«j  -  0.3)  +  0.768(w2  -  0.3)  >  6.0 

g3(x)  =  1.0  +  0.0332(x1  -  0.3)  +  0.0186  (x2  -  0.3) 
+  3.34(  x3  -  0.3)  +  0.0204IWJ  -  0.3) 
+  0.778(w2  -  0.3)  +  2.62(w3  -  0.3)  >  3.5 

0.8  <  ij  <  1.0,       i  «,  1  .  2 ,  3 
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Let  the  above  constraints  be  denoted  by  x  6  X.  Wj  is  related  to 


«j  by: 


0.39 


1.39  -  x,. 


i  »  1.2.3 


We  solve  the  problem  by  the  following  steps. 
§.iSE  2^    2klai_n  and  construct  PI.S-  &  NiS-j>ayof f  tables 
EII  -  Eayoli  Table 


1 

'l 

f2 

f3 

f4 

xl 

x2 

x3 

fl 

6 

7922 

.3413 

9.6824 

1.2110 

.9999 

.9999 

.7800 

f2 

5 

2777 

6.2752 

2.6735 

1 .9306 

.8507 

.  9090 

.8496 

f3 

4 

8599 

.3413 

1 .0418 

2 . 1468 

.9999 

.7817 

.  8638 

f4 

6 

7911 

.  3761 

9 .6824 

.9588 

.9999 

1 .0000 

.  7417 

f   =  (6.7922.  6.2752.  1.0418,  0.9588) 


NJLS    -    Payoff    Table 


fl 

4 

7843 

6 

1099 

1 

5173 

2 

0793 

.8675 

.8385 

.  8596 

f2 

5 

4931 

3413 

2 

0919 

1 

9794 

.9999 

.  8809 

.8530 

f3 

6 

4856 

5 

4252 

9 

6824 

1 

6066 

.9141 

.9999 

.  8235 

f4 

5 

3276 

5 

4192 

2 

4319 

11 

3750 

.  9144 

.  8986 

1 .0000 

f       «     (4.7843.     0.3413.     9.6824.     11.3750) 
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SteE  li    £al£Sl3ting  Phase 

We  assume  the  four  objectives  are  of  equal  importance,  i.e.  w 
«  (1.0.  1.0,  1.0.  1.0).  and  H  =  (0.5,  0.5).  We  also  assume  the 
analyst  wants   only  the  p  =  1.  2  and  oo  cases. 

The  distance  function  for  measuring  the  closeness  to  the 
positive  ideal  solution  is  given  as: 

6.7922  -  f1(x) 

F   =  Wj  [ ]   +  (3.23) 

(6.7922  -  4.7843) 

6.2725  -  f2(x)        n 
(6.2725  -  0.3413) 

f3(x)  -  1.0418         p 

w3  t 1   + 

(9.6824  -  1.0418) 

f4(x)  -  0.9588         _ 
w4  [ ]P 

(11.3750  -  0.9588) 


The  distance  function  for  measuring  the  distance  to  negative  idea 
solution . 


fjfx)  -  4.7843 
Wj  [ j 

( 6 . 7922  -  4  .  7843 ) 


(3.24) 


f,(x)  -  0.3413 

r  —  i  P 

w2  [ ] 

(6. 2725  -  0 .3413) 
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w,  [- 


9.6824  -  f,(x) 


(9 .6824  -  1 . 0418) 


w4  [- 


11.3750  -  fd(x) 


(11.3750  -  0.9588) 


Case  1  .    p  =  1 

The  problem  Is  to  find  x  so  as  to: 

min  {  Hp  Fj  -  Hn  Hj  > 

s .  t .  x  e  x 

This  is  a  nonlinear  programming  problem.  We  use  Sequential 
Unconstrainted  Minimization  Technique  (SUMT)  to  solve  the  problem. 
The  solution  is : 

x  =  (0.8507,  0.9476,  0.7726) 

f  =  (5.5624,  6.2752,  4.0285.  1.1565) 

Case  2.    p  =  2 

The  problem  is  to  find  x  so  as  to: 

min  {  Hp  F2  -  Hn  N2  > 

s  .  t  .    x  €  X 

This  is  a  nonlinear  programming  problem.  We  use  SUMT  to  solve 
it.   The  solution  is  given  as: 
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X    •     (0.8507,     0.9484.     0.7722) 

f    =     (5.6611.     6.2752.     4.0682,     1.1540) 

Case    3 .         p    =    oo 

The     problem     is     not     solvable     with     the     current     problem 
formulation.        We    rewrite     (3.23)    as: 

min    d 

s. t .       x    €    X 

6.7922    -    f , (x) 

«l     [ ]     <    d* 

6.7922    -    4.7843  (3.25) 

6.2725    -    f2(x) 

w2     [    ]     <    d* 

6.2725    -    0.3413 

f,(x)     -    1.0418 

r  O       —  * 

w3    [    ]    <    d 

9.6824    -    1.0418 

f4(x)     -    0.9588 
w4     [     ]     <    d* 

11.3750    -    0.9588 

Similarly.       we    rewrite     (3.24)    as: 

max    d" 

s  .  t .       x    €    X 


w,     [- 


fj(x)     -    4.7843 


6 . 7922    -    4 . 7843 ' 


]     >    d" 


(3.26) 
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f2  (x)  -  0.3413 
w2  [  ]  >  d 

6.2725  -  0.3413 

9.6824  -  f3(x) 
w3  [  ]  >  d~ 

9.6824  -  1.0418 

11.3750  -  f4(x) 
«4  [  ]  >  d" 

11.3750  -  0.9588 

Combination  of  (3.25)  and  (3.26)  gives  a  solvable  problem 

min  (  d*  -  d~  )  (3.27) 

s. t.    x  €  X 

6.7922  -  f i (x) 

H   Wj  [ ]  <  d* 

6.7922  -  4.7843 

6.2725  -  f2(x) 

Hp  «2  [ ]  <  d* 

6.2725  -  0.3413 

f3(x)  -  1.0418 
H   „3  [ ]  <  d 

9.6824  -  1.0418 

f4(x)  -  0.9588 

Hp  «4  t ]  <  d* 

11.3750  -  0.9588 

t1ix)     -    4.7843 
Hn  w2  [ ]  >  d" 

6.7922  -  4.7843 
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Hn  "a  (- 


f2(x)  -  0.3413 
6 . 2725  -  0.3413 

9.6824  -  f3(x) 
9.6824  -  1.0418 
11 .3750  -  f4(x) 
11 . 3750  -  0 . 9588 


>  d" 


>  <T 


]  >  d" 


The  problem  Is  solved  by  SUMT .   The  solution  is  given  as  : 
X  =  (0.9608.  0.9547.  0.8268) 
f  -  (6.0098.  3.9650.  4.4057,  1.6432) 

i.tep_  2j_      Eva_luatl,on  p_hase 

The  three  sets  of  solutions  are  summarized  in  Table  3.4.   The 
DM  may  select  one  or  none  of  them. 

Table  3.4  Solutions  Table 


f  1 


f  2 


f  3 


f  4 


x2 


PIS 

6 

7922 

6.2725 

1.0418 

0.9588 

—  - 

- 

-- 

-- 

-- 

-- 

NIS 

4 

7843 

0  .  3413 

9.6824 

11 . 3750 

- 

— 

- 

— 

- 

— 

P-l 

5 

5624 

6 . 2752 

4 .0285 

1 . 1565 

0 

8507 

0 

9476 

0 

7726 

2 

5 

6611 

6 . 2752 

4.0682 

1 . 1540 

0 

8507 

0 

9484 

0 

7722 

00 

6 

0098 

3 .9650 

4.4057 

1 .6432 

0 

9608 

0 

9547 

0 

8268 
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3  .  4   Remarks 

(1)  TOPSIS-MODM  requires  no  subjective  preference  information 
from  the  DM  once  the  problem  constraints  and  objectives  have  been 
defined.  Thus  this  method  may  require  the  DM  to  accept  the 
solution  obtained  from  the  method. 

(2)  The  solutions  obtained  in  the  calculation  phase  depend 
not  only  on  the  constraints  and  objectives  but  also  heavily  rely 
on  the  parameter  p.  Notice  that  in  our  numerical  examples,  no 
solutions  are  identical  in  any  given  solution  table.  That  is  the 
parameter  p  locates  these  solution  points  which  are  geometrically 
different  from  one  another. 

(3)  Unlike  those  methods  that  generate  a  whole  set  of 
nondominated  solutions  which  take  the  DM  a  lot  of  time  to  process 
and  absorb.  TOPSIS-MODM  generates  only  three  nondominated 
solutions.  This  feature  makes  the  selection  of  the  preferred 
solution  less  time-consuming  and  easier. 
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Chapter    4       Interactive    Sequential    Goal    Programming   II 

(ISGP    I  I    ) 
4.1        Introduction 

A  new  approach  called  Interactive  Sequential  Goal  Programming 
II  (ISGP  II)  is  introduced,  which  is  an  updated  ISGP  of  Masud  and 
Hwang  [1981].  It  combines  the  attractive  features  of  the 
following  MODM  methods:  STEM  method's  [Benayoun,  et  al..  1971] 
positive-ideal  solution  concept  and  the  interactive  calculation 
phase  and  decision  phase  process  between  the  DM  and  the  analyst; 
Charnes  and  Cooper's  [1961]  original  Goal  Programming  concept; 
SEMOPS's  [Monarchi.  et  al..  1976]  principal  problem  and  auxiliary 
problems  approach.  In  ISGP  II.  the  new  concepts  of  negative-ideal 
solutions  of  multiple  objectives  and  a  new  normalization  process 
of  each  objective  by  incorporating  the  positive-ideal  solution  and 
nega t i ve- i dea 1  solution  are  introduced.  The  Maximum  Achievable 
Rate  (MAR)  and  the  Maximum  Achievable  Goal  (MAG)  for  the  multiple 
objectives  are  also  introduced.  The  set  of  MAG  is  a  feasible 
solution    which    has    the   highest    same    rate,     MAR.     for    all    objectives. 

Maximum  Achievable  Rate  £MARi  and  Maximum  Achievable  Goai  iMAgix 
In  a  group  decision  making  problem,  if  each  objective  represents 
one  decision  maker's  payoff,  then  the  pos i t i v e- i dea  1  solution, 
f  .  if  attainable,  would  make  each  DM  happy  because  it 
simultaneously  optimizes  each  DM's  payoff.  However,  in  a  real 
world  problem,  f  is  usually  not  attainable.  To  dissolve  the 
group    conflict,    a    compromise    among   the    DM's    is    necessary.       That 
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is,  each  DM  can  expect  to  achieve  only  portion  of  the  f  instead 
of  the  entire  f  The  amount  of  sacrifice  by  each  DM  must  be 
equal  in  order  to  maintain  a  "cooperative  group  spirit".  Hence. 
MAR   and  MAG   are  introduced  to  satisfy  that  condition. 

MAG  and  MAR  are  determined  by  the  analysts  prior  to  having  the 
DM  set  initial  goals,  therefore  they  served  as  excellent  reference 
values.  The  MAG  may  be  interpreted  as  a  set  of  nondominated 
solutions  which  is  the  individual  DM's  largest  amount  of  resources 
s/he  can  have  without  the  sacrifices  of  the  other  DMs.  MAR  is  the 
maximum  percentage  of  resources  a  DM  can  have  without  jeopardizing 
the  group  harmony.  There  is  a  procedure  (see  ISGP  II  algorithm  for 
details)  which  operates  in  an  iterative  fashion  that  can  be  used 
to  derive  the  values  of  MAG  and  MAR. 

With  the  concept  of  MAR  and  MAG,  we  proceed  to  formulate  the 
algorithm  of  ISGP  II. 
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4.2    Algorithm  For  ISGP  II 

Pl2kIS5  E°E!Ed.l  at  ionx 

The  MODM  problem  is: 

Max/Min  [  f^x).   f2(x) ffclx)]  ,4'11 

s.t.   x  6  X 
where    fWx):  the  benefit  objective  for  maximization,  i  €  J 
fj(x):  the  cost  objective  for  minimization,  i  6  I 

Algorithm^ 

The  complete  ISGP  II  algorithm  for  this  problem  is  given 
below . 

Ste.E.  0    Determine  Positive  Ideal  Solution  (PIS),   f  .   and 
Negative  Ideal  Solution  (NIS).  f~. 

Steg  0-1    Obtain  and  Construct  PIS  -  Payoff  Table. 
The  problem  is  to  find  x  so  as  to 


max  fj(x)  .   j  €  J 
s.t.    x  €  X 


(4.2) 


and  /  or 


s.t.   x  6  X 


(4.3) 


The    solution    of    (4.2).     x       and    f  *     .     is    known    as    the    positive- 
ideal    solution.       Similarly,     the    solution    of    (4.3).     x      and    f ,    .     is 
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also  known  as  the  positive-ideal  solution.   With  all  the  x   from 
(4.2)  and  (4.3).  the  PIS-Payoff  Table  is  given  as: 

Eli  -   Eay°ll  !akl£ 


L21 


12 


■  Ik 


2k 


lkl 


Lk2 


f   =  (  f 


) 


Step.  0-2   Obtain  and  Construct  NIS-Payoff  Table. 
The  problem  is  to  find  x  so  as  to 


min  f  •  ( x )  .   j  €  J 
s. t.    x  €  X 


(4.4) 


and  /  or 


max  f  j (x) .  i  6  I 
s. t.    x  €  X 


(4.5) 


The  solutions  of  (4.4)  and  (4.5).  x~  and  fr~-  are  known  as  the 
nega  t  i  ve-  idea  1  solutions.  With  all  the  x"  from  (4.4)  and  (4.5). 
the  NIS-Payoff  Table  is  constructed  as: 


39 


NIS-Payoff  Table 


1  ___2  _ 
1  z12 
21     f  2 


zlk 
z2k 


akl 


2k2 


f 


) 


Note  that  f  and  f~  provide  the  upper  and  lower  bounds  on  the 
objective  function  values  and  are  a  great  source  of  information 
for  the  DM  in  all  the  following  steps. 

Step  1   Set  initial  goals. 

Step.  1.1   Obtain  a  set  of  Maximum  Achievable  Goals  (MAG)  and  a 
Maximum  Achievable  Rate  (MAR).  B%.     (The  MAG  is  a  feasible  solution 
which  has  the  highest  same  rate,  MAR.  for  all  objectives.) 
The  problem  is  to  find  x  so  as  to 


{  £   W_  d.   > 

r=i   r   r 


(4.6) 


s. t.    x  6  X 


fj(x)  +  (dj 


f  .  (x)  -  (d f 


dJ  »  (  fj 


(  f. 


fj"  > 


J  €  J 


i  6  I 


(4.6a) 


(4.6b) 


0  <  d. 


1.2. 
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where 


bi 


fj        -     (1     -    B*)     (fj 
f4*    +    (1     -     B*)    (f  ,  " 


fj    ),     j    €   J 


f,     )  .    16    I 


The    solution    is    given    as    x.     fg.     and    ag. 

Note    that     (4.6a)    and    (4.6b)    are.     respectively,     transformed 
from 


and 


fjfx) 


-    f. 


fill) 


(d. 


dJ  » 


(d. 


j 


j    6   J 


i    €    I 


Also  note  that  MAR  is  unknown  initially.  To  properly  start  the 
calculation  of  (4.6),  we  have  to  assign  a  value  to  B%.  Since  we 
perceive  that  when  B*  is  set  at  75*  we  can  obtain  a  MAR  with 
minimum  effort,  we  suggest  B*  be  set  at  75*  to  start  the 
calculation      of      (4.6).  Steps      1.1.1      and      1.1.2      illustrate      the 

detai  Is . 

Step_    1.1.1       Check    if    fg    is    a    set    of    nondominated    solutions. 
If    ag    -    0.     go    to    step    1.1.2. 

If  ag  ^0,  we  decrease  B*  by  5*.  i.e.  B%  =  (B  -  5)*, 
recalculate  (4.6).  and  check  the  value  of  ag  again.  The  procedure 
is  repeated  until  an  answer  ag  =  0  is  found.  After  an  as  =  0  is 
found,     we    increase       BX    by    1*.     i.e.     B*    =    (B    +     1)*.  recalculate 

(4.6).     and    check    the    value    of    ag    again.       The    procedure    repeats 

/ 
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until  an  ag  4  0  is  found,  and  that  a^  4  0  signals  that  MAG  and  MAR 
are  found.  MAR  is  the  highest  rate  for  all  objectives  in  (4.6)  for 
which  a  feasible  nondomi nated  solution  is  obtainable;  this  is  MAG. 
Go    to    step    1.2     . 

Step_    Ij_L.£ 

When  ag  =  0.  we  increase  B%  by  5*.  i.e.  BX  =  (B  +  5)%. 
recalculate  (4.6),  and  check  the  value  of  a*>  again.  The  procedure 
is  repeated  until  an  answer  ag  4  0  is  found.  After  an  ag  4  0  is 
obtained,  we  decrease  B9>  by  1*,  i.e.  B*  -  (B  -  l)X.  recalculated 
(4.6).  and  check  the  value  of  ag  again.  The  procedure  is  repeated 
until  an  answer  ag  =  0  is  found,  and  that  ag  =  0  signals  that  MAG 
and  MAR  is  found.  We  conclude  that  MAR  is  (B  +  1)%  and  MAG  is  a 
set  of  nondom i nated  solutions  obtained  at  (B  +  1)*.  Go  to  Step 
1.2.  Notice  that  MAG  and  MAR  are  used  as  another  important 
reference  in  all  the  following  steps. 

Step.  1^2    Display  ISGP  II-Initial  Table. 


max/min 


max/min 
f2(x) 


max/min 


f 

MAG  (B*) 
f 

w 


The  initial  table  is  presented  to  the  DM.    If  the  DM  is 
satisfied  with  the  MAG  then  take  it  as   the  preferred  solution  and 
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go  to  step  8.   Otherwise  go  to  step  1.3  . 

Step.  l_j_3    Set  initial  goals. 

The  DM  set  his/her  initial  goals  at  b  referring  to  the  initial 

table.  The  initial  goals,  b.  should  satisfy 


fj"  <  °i  -  fi  '  J  €  J 


(4.7) 


fi   <  bj  <  fj  ,  i  €  I 

Step.    2       Solve    for    PS    (Principal    Solution). 

The     goal     constraints     associated     with     the     objectives  are: 


and 


fj       -    fj(x) 


fj       -    f  J 


(dj       -    dj     )    - 


f J       -    bJ 


fJ       -    f i 


1    €    J 


(4.8) 


(d.       -    d,     ) 


i    6    I 


such    that    (4.7)    will    not    be    violated. 

Expression     (4.8)     can    be     rewritten    as: 


fj(x)     +     (dj       -    djT)    (fj 


fj~' 


J    6    J 


di+>    (fi"    "    fi    >     "    bi-     i    e    l 


(4.9) 
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££SE  i_J.   Calculation  of  PS  (Principal  Solution). 
The  problem  is  to  find  a  set  of  x  so  as  to 

min  aps  =  {  £  Wr  dr~  }  (4.10) 

s. t.    x  6  X 

fj(x)  +  (dj-  -  dj+)  (fj*  -  fj")  -  bj.  J  €  J 
f^x)  -  (dj-  -  di+)  ltx-    -  fj*)  -  bA.  16  1 


0  <  dr  ,  dr+  <  1 


where  W   is  the  relative  importance  of  each  objective  assigned  by 

the  DM. 

Solution:   x° .  f°.  a„0 
—    —     p  s 

d  ~  is  the  underachie vement  of  the  objective  function  f_(x)  from 
its  goal  b_.  £   W   d  ~  equals  zero  if  the  goal  is  in  the   feasible 


the  feasible  region. 

Steg  2.2  Check  if  f°  is  a  set  of  nondominated  solutions 

If  a„„  «  0  then  go  to  step  2.3. 
v  s 

If  a  4    0  then  take  f°  as  the  principal  solution  and  go  to 
step  3 . 

Ste£  2_;_3  Improvement  of  f  ° . 


PS 


-  0  implies  that  the  f°  obtained  is  a  set  of  dominated 


solutions.   We  increase  bj  and  decrease  bj  by  z% ,     recalculate 
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(4.10),  and  check  the  value  of  a 


ps  • 


The  procedure  repeats  until 


an  a  /  0  is  found.  We  recommed  z  to  be  set  at  5.  10.  15.  or  20* 
of  the  goal  of  each  objective  with  a  default  of  5*.  Notice  that 
the  goal  setting  rules  in  (4.7)  must  be  satisfied  when  modifying 
b  •  and  b  j . 

§.t£P_  2    Solve  for  k  sets  of  ASs  (Auxiliary  Solutions). 

Each  set  of  auxiliary  solutions  is  obtained  by  setting  one 
individual  goal  as  a  constraint  and  by  minimizing  the  summarized 
weighted  deviations  from  other  goals. 

The  problem  is  to  find  a  set  of  x  so  as  to 


min  au  =  {  Z   Wr  d~  } 

UK 

s.t.    x  €  X 


(4.11) 


fj(x)  +  (dj~  -  d^)  (fj   -  fj  )  =  bj,  (for  all  j.  j  4    u) 
fjlx)  -  (dj"  -  di*)     (fj"  -  fi*)  =  bj,  (for  all  i,  i  4    u) 


f  u(£>  1  bu   (if  u  €  I) 

0  <  dr~  <  1     (for  all  r,  r  t    u) 

where  Wr  is  the  relative  importance  of  each  objective  assigned  by 

the  DM. 

Solution:    xu.  fu.  au.  u  =  1.2 k 
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§i£E  3- 1    Check  if  fu  is  a  set  of  nondominated  solutions. 

If  a"  =  0  then  go  to  step  3.2. 

If  au  ±  0.  then  we  take  fu  as  the  utn  auxiliary  solution  and 
go  back  to  step  3  for  the  (u  +  1  )  auxiliary  solution.  If  all  k 
sets  of  auxiliary  solutions  are  obtained,  i.e.  u  +  1  >  k,  we 
proceed  to  step  4. 

Step_  3.2    Improvement  of  fu. 

au  =  0  implies  that  the  previously  obtained  fu  is  a  set  of 
dominated  solutions.  We  increase  bj  and  decrease  b^  by  z% , 
recalculate  (4.11).  and  check  the  value  of  au  again.  The 
procedure  repeats  until  an  au  4  0  is  found,  i.e.  a  set  of 
nondominated  solutions  is  obtained.  Z%  has  the  same  featuresas 
those  in  step  2.3  .  Notice  that  the  goal  setting  rules  in  (4.7) 
must  be  satisfied  when  modifying  b*  and  b(.   Go  to  step  3. 


Step.  4    Form  ISGP  II-Tradeoff  Table. 

The  tradeoff  table  is  given  in  Table  4.1.  The  first  row 
marked  f  represents  the  positive-ideal  solutions.  The  second  row 
marked  MAG  (B9>)  indicates  the  maximum  achievable  goal  and  the 
maximum  achievable  rate.  The  third  row  marked  f~  represents  the 
negative  ideal.  The  fourth  row  marked  w  indicates  the  weights  of 
each  objective  assigned  by  the  DM.  All  four  rows  appear  on  the 
tradeoff  table  for  each  iteration.  The  notation  "#"  right  below 
dotted  line  indicates  the  number  of  iterations.  The  fifth  row 
marked  Goal  indicates  the  goals  set  by  the  DM.   The  Sixth  row 

/ 

96 


marked  PS  Indicates  the  current  best  nondominated  solution, 
of  the  remaining  rows  indicates  a  different  AS. 


Each 


Table  4.1   ISGP  II-Tradeoff  Table 


max/min       max/min      ...     max/min 
fjlx)         f2<x)               fk(^> 

*         *             *                 „  • 
f        fj            f2          ...      fk 

MAG  (B*)      bj*           b2*          ...      bk 

f"       tf                            f2"          ...       fk" 

w        Wj            «2                         •  •  •      wk 

Iteration  # 

goal         bj            b2          .  .  .      bk 

PS         fj°          f2°        ...      fk° 

ASl          bx            f2l                       ...      fk: 

AS2          fj2           b2          ...      fk2 

ASK 


§i£E  5    Evaluation  phase. 

Step.  5^1   Ask  the  DM: 
Is  f°  or  any  fu  satisfactory  ? 
If  yes.  then  go  to  step  8.  otherwise  go  to  step  5.2. 

Stej>  5.2  Ask  the  DM  to  specify  his/her  new  goals,  b.   from  the 
Information  in  the  tradeoff  table. 
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Step.  6   Calculate  the  deviations  of  the  new  goals  from  PS. 
tf    =  fj°  -  bj,  for  all  J  =  {  j  |  bj  <  fj°  > 


'j  -  bJ 


for  all  j  -  {  j  |  bj  >  fj°  ) 


?i~  "  bi  "  fi°'  for  all  1  =  {  i  I  bj  >  fi0  > 
?i  +  -  fi°  "  bi'  for  ^1  1  =  {  i  !  bj  <  fj0  } 

St eg  7   Consistency  check  for  the  specified  b. 

For  all  r  =  1.2 k,  check  to  see  if  there  is  at  least  one 

A   _  A   + 

f    >  0  when  one  or  more   fr   >  0. 

If  yes.  then  go  to  step  2.1.  otherwise  go  to  step  5.2  .  Since 
f°  is  a  nondominated  solution,  this  consistency  check  is  necessary 
to  make  sure  that  if  the  DM  wants  an  improvement  in  one  of  the 
objective  values  from  f°.  s/he  is  ready  to  accept  a  detriment  in 
at  least  one  other  objective  value. 

St £2  8    Stop. 

Print  the  preferred  solution  x  and  f  and  terminate  the 
algorithm . 

Note  that  in  solving  (4.6),  (4.10).  and  (4.11),  any  appropriate 
single  objective  optimization  technique  can  be  utilized 
i teratively . 
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4.3   Numerical  Examples 

4.3.1     Hardee  Company  Production  Scheduling. 

The  Hardee  company  makes  two  kinds  of  dolls.  Doll  A  is  a  high 
quality  doll  and  Doll  B  is  of  lower  quality.  The  respective 
profits  are  $0.40  and  $0.30  per  doll.  Each  Doll  A  requires  twice 
as  much  labor  as  a  Doll  B.  and  if  all  dolls  were  of  type  B,  the 
company  could  make  500  per  day.  The  supply  of  material  is 
sufficient  for  only  400  dolls  per  day  (both  A  and  B  combined). 
The  problem  assumes  that  all  the  dolls  for  type  A  and  type  B  the 
factory  can  make  could  be  sold,  and  that  the  best  customer  of  the 
company  wishes  to  have  as  many  as  possible  of  type  A  doll.  The 
manager  realizes  that  two  objectives:  (1)  the  maximization  of 
profit,  and  (2)  the  maximization  of  production  of  Doll  A.  should 
be  considered  in  scheduling  the  production. 

This  M0DM  problem  can  be  formulated  as  follows: 


max  f  ,  (  x  )  =  0.4x 


max  f , ( x ) 


1  + 


0  .3x. 


subject  to 


,(x)  »  2  Xj  +  x2  <  500 

Xj .  Xj  >  0 


where  Xj  and  x2   are  the  number  of  Doll  A  and  that  of  Doll 
produced  respectively.   Let  constraints  be  denoted  by   x  €  X. 
The  problem  is  solved  by  the  following  steps. 
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Step.  0   Obtain  and  Construct  PIS  &  NIS  Payoff  Tables. 


We  solve  max  fjfx)  and  max  f2'x)  subject  to  the  constraints, 
respectively,  and  obtain  the  PIS-Payoff  Table  as: 


130 
100 


100 

* 
250 


f 


(  130 .  250) 


Similarly,  we  solve  min  fj(x)   and   rain  f2(x)  subject  to   x   X. 
respectively,  and  obtain  the   NIS-Payoff  Table  as: 


t~       =     (0,     0) 

Step.    UL         Obtain    a    set    of    MAG    and    MAR. 

We    (the    analysts)    set       B    *    =    75*.     and    solve    the    following 
problem : 


min    ag    =    (    dj       +    d2       ) 


s.t.      x     €    X 


(4.5) 


100 


f j (x)  +  130  (d1~    -     d1+)    -  bj 

f  2(x)  +  250  (d2~  -  d2+)  =  b2 

0  <  dr"  .  dr+  <  1  .  r  =  1.  2 

where  bj  =    130  -  (  1  -  75*  )  (  130  -  0  )  =  97.50 
b2  •  250  -  (  1  -  75%  )  (  250  -  0  )  -  187.50 

This  is  a  linear  programming  problem.   The  solution  is  given  as: 
x  =  (187.50.  75.00) 
f  g  =  (97.50.  187.50) 

ag  =  0 
The  answer  ag  =  0  implies  that  with  the  given  b  there  exists  a 
dominated  fg   which  is  not  appropriate  to  serve  as  MAG. 

Since  B%  ■  75%  results  in  a  dominated  b   and  fg  .  we  increase 
B%   by  5%.    i.e.  B%   =  80%.  and  recalculate  (4.5)  with 

bj   =  130  -  (  1  -  80%  )  (  130  -  0  )  =  104.00 
b2   ■  250  -  (  1  -  80%  )  (  250  -  0  )  =  200.00 

The  solution  is  given  as 

x  =  ( 200.00  .  80.00) 

fg  =  (104.00.  200.00) 

ag    =    0 
Again,      we     find     the     new    ag     =     0.      i.e.      fg     is     a     set     of     dominated 
solutions.        We     can     increase     B%    by     another     5%.      i.e.     B%     ■     85%. 
Another    calculation    of    (4.5)    is    conducted    with 
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35*    )    (130    -     0)    =    110.50 
(     1     -     85*    )     (250    -     0)    =    212.50 

The    solution    is    given   as: 

x    =    (212.50.     75.00) 
f g    =    (107.50.    212.50) 
ag    -    0.0231 

The  answer  ag  t  0  implies  that  fg  is  a  nondominated  solution.  The 
MAR,  B*,  is  located  somewhere  between  85%  and  81*.  To  locate  the 
precise  B*.  we  decrease  B*  by  1*.  i.e.  B*  ■  84*.  and  calculate 
(4.5)  with  bx  =  109.20  ar 
The  solution  is  given  as 

x  ■  (210.00,  80.00) 

f  g  •  (108.00,  210.00) 

ag  =  0.0092 
Since  the  answer  ag  is  not  zero,  we  decrease  B%  by  another  1*. 
i.e.  B*  -  83*.  and  proceed  to  another  calculation  of  (4.5)  with   b 
=  (107.90.  207.50).   The  solution  is  given  as 
x  =  (207.50.  83.00) 
f  g  =  (107.90.  207.50) 
ag  =  0 
The  answer  ag  -  0   Implies  that  the   fg  when  b  ■  83*  is  a  set  of 
dominated  solutions.   We  conclude  that  MAG  is  located  at  fg  ■ 
(108.00,  210.00)  and  MAR  Is  84*.  i.e.  B*  -  84*. 
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The  scheme  for  obtaining  MAG  and  MAR  is  depicted  in  figure 
4.1.  Points  P  and  N  denote  the  positive-ideal  solution  (130.00. 
250.00)  and  nega t i v e - i dea 1  solution  (0.  0),  respectively.  The 
heavy  solid  line  AB  represents  the  nondominated  solution  region. 
Capital  S  denotes  the  feasible  region  for  our  problem.  Point  1.  f 
=  (97.50.  187.50).  is  the  set  of  solutions  when  B*  =  75*  (or  b  - 
97.50,  187.50)).  Point  2,  f  =  (104.00.  200.00).  is  the  set  of 
solutions  when  B*  =  80*  (or  b  =  (104.00.  200.00)).  Point  3.  f  = 
(107.50.  212.50).  is  the  set  of  solutions  when  B*  =  85*  (or  b  • 
(110.50.  212.50)).  Points  4  and  5  represent  sets  of  solutions 
when  B*  =  84*  and  B*  =  83*.  respectively.  It  is  obvious  that 
points  1.2,  and  5  are  dominated  solutions  while  points  3  and  4  are 
nondominated  solutions  (they  are  on  the  AB  line).  Since  point  4 
is  derived  from  a  more  precise  B*.  B*  -  84*.  than  that  of  B*  = 
85*.  we  choose  point  4.  f  =  (108.00.  210.00).  as  MAG  and  84*  as 
MAR  for  our  problem. 

§.tep_  1^2    Display  ISGP  II  -  Initial  Table. 


Initial  Table 


max  fjfx) 


max 

'l' 

[£) 

250 

00 

210 

.00 

0 

.00 

1 

.00 

MAGI84*) 
f" 

w 


130.00 

108.00 

0.00 

1  .00 


The  initial  table  is  presented  to  the  DM.   We  assume  the  DM  does 
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Fig.  4.1  The  procurement  of  MAR  and  MAG  for  Hardee 
Co.  production  schedule  problem. 
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not  accept  the  MAG  (84*)  as  the  preferred  solution. 

Step.  l._J3     Set  Initial  Goals  b. 

Assume   the  DM   sets   his/her   initial   goals   b   at   90%    of 
fj*  -  f  1~     I.  i.e.  b  =  (117.00.  225.00).  referring  to   the  ISGP 
II  initial  tabl  e. 

lists*  i.22  Nck_  1^ 

Stg_2  2A.        Solve  for  PS  (Principal  Solution). 
The  problem  is  to  find  x  so  as  to 
min  aps  =  (dj"  +  d2"> 

s. t.    x  €  X 


fj(x)  +  130  (dj   -  dj+)  =  ba  =  117 


J2   - 

This    is    a    linear    programming    problem.       The    solution    is    given    as 
X    =     (225.00,     50.00) 
| °    =    (105.00.     225.00) 


Ste£    2^.2 

As    aps    ^     °    imPlies    that    the    solution    f°    is    a    nondominated 
solution    and    is    accepted    as    the    principal    solution.       Go    to    step    3. 
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Step  3    Solve  for  ASs  (Auxiliary  Solutions). 

The  first  auxiliary  solution  is  obtained  by  solving 

min   a   =  {  d?   } 


fj  (x)  >  bj  =  117 

f2(x)  +  250  (d2"  -  d2+)  =  b2  =  225 


0  <  d2  .  d2+   <  1 


This  is  a  linear  programming  problem.   The  solution  is  given  as 
x  =  (165 .00.  170.00) 
f1  =  (117.00.  165.00) 
a1  jt    0 

Step.  3^1 
The  ans 

accepted  as  the  first  auxiliary  solution.   Go  to  step  3  for  the 
second  auxiliary  solution  set. 

Step  3   Next,  we  obtain  the  second  auxiliary  solution  by  solving 

min  a2  =  (  dj"  } 

s.t.   x  £  X 


f1(x)  +  130(dj   -  A 1  * )     =bj=  117 
f,  (x)  >  b,  =  225 
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0  <  dj  .  dj+  <  1 


This  is  a  LP  problem.   The  solution  is  given  as 
X  =  (225.00  .  50.00) 
f  2  =  (105.00.  225.00) 
a2  f    0 

Ste^  3. 1 

Since  a2  4    0  we  accept  f   as  the  second  auxiliary  solution. 


Step  4    Form  ISGP  II  -  Tradeoff  Tablf 
max  f,  (x) 


max  f 2 (x ) 
250.00 
210 .00 

0  .00 

1  .00 


f* 

130 

.  00 

MAG (84*) 

108 

.  00 

f~ 

0 

.  00 

w 

1. 

.00 

Iteration  1 

Goal 

117 

.  00 

PS 

105 

00 

AS1 

117 

-  00 

AS2 

105 

,  00 

225 .00 
225.00 
165.00 
225.00 


This  tradeoff  table  is  presented  to  the  DM  for  evaluation. 


Step.  5     Evaluation  phase. 

We  assume  the  DM  is  satisfied  with  none  of  the  solutions, 
proceed  to  Step  5.2  for  setting  another  set  of  goals  b. 


We 
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Step.  5^2    Set  a  set  of  new  goals  b. 

Since  none  of  the  PS  and  ASs  are  satisfactory  to  the  DM.  a  set 
of  new  goals  b  is  set  at  (110.00.  215.00)  by  the  DM. 

*  -        A  + 
Ste£  6    Calculation  of   f     and   f   . 

The  new  goals,  b  ■  (110.00.  215.00)  are  compared  with  f°  « 

(105.00,  225.00)  of  the  first  iteration. 

A   4- 

fl 

A   _ 

f2 

Ste£  7    Consistency  check  on  new  b. 

There  is  one  f  ~  >  0   and   one  f+  >  0.   That  is.  the  new 
goals  b  is  a  set  satisfying  the  consistent  test  and  is  accepted  as 
the  initial  goals  for  the  second  iteration.   Go  back  to  Step  2.1 

I£S.E£t  j.on  No. 2 

SteE    <L_1         Solve    for    PS    (Principal    Solution). 

The   problem   of    (4.6)    is   solved   again   with   bj    =    110.00    and      b2 
"    215.00    .       The    solution    is    given    as 
X    •     (215.00,     70.00) 
f °    -     ( 107.00 .     2  15.00) 
a"s    fl    0 

Step.    2^         The    answer    aps     4     0       implies     that     f°     is     a    set    of 
nondominated    solutions    and    is    accepted    as    the    principal     solution. 
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Go  to  Step  3. 

Steg  3    Solve  for  ASs  (Auxiliary  Solutions). 

Similar  to  that  of  iteration  1,   We  solve  (4.7)  with  bj  » 
110.00  and  b2  ■  215.00.    The  solution  is  given  as 

x  =■  (200.00  .  100.00) 

f1  =  (110.00.  200.00) 

a1    jt    0 

Step_  3 . 1  The  answer  a1  4  0  implies  that  f1  is  a  set  of 
nondomi nat ed  solutions  and  Is  accepted  as  the  first  auxiliary 
solution.   Go  to  Step  3. 

S£<LE  3    This  step  solves  for  the  second  auxiliary  solution. 

Problem  (4.8)  was  s 
The  solution  is  given  as 

x  =  (215.00,  70.00) 

f  2  =  (107.00.  215.00) 

a2  t    0 

Si£E  3^1  The  answer  a2  ^0  implies  that  f2  is  a  set  of 
nondominated  solutions  and  is  accepted  as  the  second  auxiliary 
solution.  Since  both  goals  have  been  satisfied  serially,  we 
proceed    to    Step    4. 

Step_    4         Form    the    ISGP    II    -    Tradeoff    Table. 
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ISGP  II  -  Tradeoff  Table 


max  fj(x) 
130.00 
108.00 

0.00 

1.00 


max  f  2 ( x ) 

250.00 

210.00 

0.00 

1.00 


f 

MAGI84*) 

f  " 

w 


Iteration  2 

Goal 

PS 

AS1 

AS2 


110.00 
107.00 
110.00 
107.00 


215.00 
215.00 
200.00 
215.00 


This  tradeoff  table  is  presented  to  the  DM  for  evaluation. 

Step.  5    Evaluation  phase. 

We  assume  the  DM  is  satisfied  with  the  first  auxiliary  solution 
AS1  and  selects  it  as  the  preferred  solution.   Go  to  Step  8. 

St  eg  8   Stop. 

The  final  solution  is  printed  as 

x  =  ( 200.00 .  100.00 ) 

f    =    (110.00 ,     200.00) 
and    the    algorithm    terminates    here. 
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4.3.2    Nutrition  Problem. 

This  problem  is  to  determine  the  quantities  of  six  foods  that 
should  be  eaten  to  balance  nutritional  requirements  so  as  to 
(i)  maximize  carbohydrate  intake, 
(ii)  minimize  cholesterol  intake,  and 
(iii)  minimize  cost. 

Information  on  the  six  foods  is  given  in  the  following  table. 
lakl£  ij.1  Nutrition  and  cost  information  of  six  foods 

Lettuce   Orange   daily 

Milk   Beef   Eggs  Bread     salad     juice   intake 

for 
( pint )( pound) ( dozen )( ounce )( ounce  )  (pint)   adults 


Vi  tamin  A  ( i . u . ) 

720 

107 

7080 

0 

134 

1000 

5000 

Food  energy 

344 

1460 

1040 

75 

17 

.4 

240 

2500 

( calories  ) 

Cholesterol 

10 

20 

120 

0 

0 

0 

(unit) 

Protein  (g) 

18 

151 

78 

2. 

5 

0 

2 

4 

6  3 

Carbohydrate ( g) 

24 

27 

0 

15 

l . 

1 

52 

Iron  (mg) 

0. 

2   10  . 

1   13. 

2 

0  . 

75 

0 

15 

1  . 

2 

12.5 

Cost  (S) 

0  . 

22   2  . 

2    0  . 

8 

0  . 

1 

0. 

05 

0. 

26 

The  problem  is  summarized  as 

max  f j (x )  =  24x j +  2 7 x 2 *  1  5x4 +  1 . 1 x g +  5 2 x 6  (  Car bohy dra t e  intake) 
min  f2(x)  ■  10xj+  20x2  +  120x3  (Cholesterol  intake) 

min  f3(x)  ■  0.22XJ+  2.2x2+  O.8X3+  0.1x4+  0 . 05x5+  0.26xg  (Cost) 
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Subject  to: 

720xj  +  107x2  *  7080x3  +  134x5  +  1000x6  >  5000   (Vitamin  A.i.u.) 
0.2XJ  +  l°-lx2  +  l3-2x3  +  0.75x4  +  0.15x5+  1 . 2xg  >12.5  (Iron.mg) 

t5+  240x6 

<6 


x.  <  6.0 


X,  <  0.25 


xR  <  4.0 


(Milk,  pints) 
( Bef f , pounds  I 
(Eggs ,  dozen) 
(Bread ,  ounces  ) 
(Lettuce  &  Salad  .  ounces  ) 
(Orange  juice,  pints) 


J  =  <1>.   I  =  {2.  3} 
The  problem  is  solved  by  the  following  steps. 
SteE  0    Obtain  and  Construct  PIS  &    NIS  Payoff  Tables. 
PIS-Payoff  Table 


1 


f , 


'1 


max  f  j 

540 

* 
00 

80.00 

6  .09 

6.  00 

1  .00 

0.00 

10 

00 

10.  00 

4 

00 

min  f  2 

380 

39 

8.44* 

3.47 

0  .00 

0  .42 

0  .  00 

10 

00 

10  .00 

4 

00 

min  f Q 

281 

67 

67.83 

2  .  26* 

3  .  78 

0.  00 

0.25 

10 

00 

0  .  00 

0 

7  9 

f   =  (540 .00 .  8.44.  2 . 26) 
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NIS-  Payoff  Table 


max  f:   93.34    71.53   3.08 


min  f,  321 .00  110. 00    4  .  75 


mln  f,  540.00  110.00   6.29" 


2  .50  0 .83  0 . 25  0. 00 
6 .00  1 .00  0.25  10 .00 
6 .00   1 .00   0 . 25   10.00 


10 .00   0 .00 

0 . 00   0.00 

10 .00   4 . 00 


f   =  (93. 34 .  110.00.  6.29) 

Step.  1^1  Obtain  a  set  of   MAG  and  MAR. 
Set   B*  •    75*  (initially) 
The  problem  is  to  find  a  set  of  x  so  as  to 


(dj  +  d2   +  d3   } 


s  .  t .  x  6  X 


fj(x)  +  (dj"  -  d/)  (fa   -  fj  )  "  bj 


f,(x)  -  (d,   -  d9+)  (f2   -  f,   )  =  b2 


f,(x)  -  (d  "  -  d,+)  (f,"  -  f,  )  -  b< 


0  <  dr".  dr+  <  1   .  r  =  1.  2,  3 


where  bj  «  540.00  -  (  1  -  75*  )  (  540.00  -  93.34  )  «  428.33 
b2  =  8.44  +  (  1  -  75*  )  (  110.00  -  8.44  )  =  33.83 
b3  =  2.26  +  (  1  -  75*  )  (  6.29  -  2.26  )  =  3.27 

This  is  a  LP  problem.   The  solution  is  given  as 
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x  =  (2.8060.  0.2885.  0.  10.  0.  3.9077) 
fg  =  (428.33.  33.83.  3.27) 
ag  =  0 

Since  ag  =0.  fg  =  (428.33,  33.83,  3.27)  is  a  set  of  dominated 
solutions  which  is  not  qualified  to  serve  as  the  problem's  MAG. 
Therefore,  we  increase  the  B*  by  5*.  i.e.  B*  «  80*.  and  repeat  the 
calculation  of  (4.9)  and  (4.10)  and  obtain  the  following  answers: 

bj  =  450.67.  b2  •  28.75.  b3  «  3.07 
X  -  (2.8752.  0.  0.  10.00.  7.8691.  4.00) 
ts    =     (435 .66,  28. 75,  3.07) 
ag  =  0.336 

Since  ag  4  0,  fg  -  (435.66.  28.75,  3.07)  is  a  nondominated 
solution.  The  MAR.  B*.  is  in  the  range  between  76*  and  80*.  To 
identify  the  precise  B*  value,  the  following  calculation  is 
conducted . 

We  decreaseB«  by  1*.  i.e.    B*  =  79*.   calculate  (4.9)  and 
(4.10),  and  obtain  the  following  solutions: 

bx  =  446.20,  b2  =  29.77,  b3  -  3.11 

X  =  (2.9768.  0.  0.  10.00.  8.2221.  4.00) 

fg  -  (438.49,  29 . 77  .  3.11) 

ag  =  0. 0173 

Since  ag  4  0.  we  decrease  B%  by  1*  again,  i.e.  B*  =  78*.  and  solve 
(4.9)  and  (4.10).  The  solutions  are  : 

/ 
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3.15 

X    •     (3.0780,     0.     0.     10.00.     8.6568.     4.00) 
Is    =    (44  1.39,     30.78.     3.15) 
aB   =    0.00096 

Since    ag    4    0.    we    reduce    the      B*   from    78*   to    77*.    and    solve    (4.9) 
and    (4.10)    again.       The    solutions    are 


'1 


437.27.     b2    =    31.80.     t>3    =    3.1! 


X    =    (2.9491,     0.1154.     0.     10.00.     4.8857.     4.00) 
f s    -    (437.27 ,     31.80.     3.19) 
ag    =    0.00 

Since  a8  =  0.  we  conclude  that  the  MAR  is  78*.  i.e.   B*  -  78*.  and 
the  MAG  is  (441.39.  30.78.  3.15). 

Steg  lj.1   Display  ISGP  II  -  Initial  Table. 


max  f , 


min  f a 


m  i  n  f  -. 


f  540.00  8.44 

MAGI78*)  441.39  30.78 

f"  93.34  110.00 

w  1 . 00  1 .00 


2  .  26 

3  .  15 
6  .  29 
1  .00 


This  initial  table  is  presented  to  the  DM.  If  the  DM  is 
satisfied  with  the  MAR  (78*),  then  take  the  MAG  as  the  preferred 
solution  and  go  to  step  8,  otherwise  go  to  step  1.3  . 
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I£££2tion  M2^il 

Step.  1^3   Set  initial  goals  b. 

We  assume  that  the  DM  set  his/her  initial  b  at  90*  of  |  fj 
fi~|.  i-e.   b»  (495.33.  18.60.  2.66).  referring  to  initial  table. 

Step.  2.1   Calculate  Principal  Solution. 

The  problem  is  to  find  a  set  of  x  so  as  to 


■in  aps  =  {dj" 


d3-  } 


x  €  X 
fj(x)  +  (d^-  d1  +  )     (fi*-  fl~)  =  bl  =495.33 
f2(x)  -  (d2"  -  d2+)  (f2"  -  fa*)  ■  b2  =  18.60 
f3(x)  -  (d3"  -  d3+)  (f3~  -  f3*)  "  b3  =  2.66 
d3-  <  1 

This  is  a  LP  problem.   The  solution  is  : 

x  «  (2.1433.  0.  0.  10.00.  3.0196.  4.00) 

f°  =  (412.77.  21.43.  2 .663) 

aps  •  0.0225 
The  answer  anc,¥     0   indicates  that  f°  is  a  set  of  nondo  m  i  na  t  ed 

pa  — 

solutions;   we  take  f  °  as  the  principal  solution. 

Step,  3   Calculate  the  first  Auxiliary  Solution  (AS1). 
The  problem  is  to  find  a  set  of  x  so  as  to 
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min  a   =  {  d2  *    ^3       ) 


s  .  t .    x  6  X 

fj(x)  >  ba  -  495.33 

f2(x)  -  (d2"  -  d2+)  (f2_  -  f2*>  •  b2  =  18.60 

f3(x)  -  <d3-  -  d3+)  (f3-  -  f3*)  -  b3  =  2.66 

0  <  d2~.  d3_  <  1 

This  is  a  LP  problem.    The  solution  is 
x  -  (5.7222.  0,  0,  10.00.  0.  4.00) 
f1  •  (495.33.  57.22.  3.30) 
a1  -  0.5381 


Step_  3^1 

a1  4    0  indicates  that  t}    is  a  set  of  nondomina ted  solutions: 
we  take  f1  as  the  first  auxiliary  solution.   Go  to  step  3. 

Step_  3   Calculate  the  second  Auxiliary  Solution  (AS2). 
The  problem  is  to  find  x  so  as  to 


Min  a2  =  <  dj   +  d3   ) 


t  .   x  6  X 


f^x)  ♦  (dj"  -  dj+)  (fj   "  tX     >  "  bl  '  495.33 


f,(x)  -  (d,~  -  d,+  )  (f,   -  f,  )  -  b3  =  2.66 


f2 (x)  <  b2  »  18 .60 


117 


0  <  d:  .  d3   <  1 


This  is  a  LP  problem.    The  solution  is  : 

X  -  (1.8596.  0.  0,  10.00.  8.6378.  4.00) 
f2  =  (412.13.  18.60.  2.88) 
a2  -  0.2404 


Step.  3^.1 

a   t*  0  indicates  that  f   is  a  set  of  nondominated  solution;  we 
take  f2  as  the  second  auxiliary  solution.   Go  to  step  3 

Step  2   Calculate  the  third  Auxiliary  Solution  (AS3). 
The  problem  is  to  find  a  set  of  x  so  as  to 


{  d.   +  d,   } 


s.t.    x  €  X 


f  !<£)+  (df-  d1*)(t  j  -  f  j  )=  bj  = 


=  4  9  5.33 


f2(x)  -  (d2"  -  d2+)  (f2_  -  t2    )   '   bj  -  18.60 


f 3 (x)  <  b3  =  2.66 


0  <  dj  .  d2   <  1 


This  is  a  LP  problem.    The  solution  is: 

X  -  (2.1433.  0.  0.  10.00.  3.0296,  4.00) 
f 3  =  (412 . 77 .  21 . 43 .  2 .66) 
a3  =  0 . 2128 
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Step.  3^1 


4    0  indicates  that  f°  is  a  set  of  nondomi nat ed  solutions; 


we  take  fJ  as  the  third  auxiliary  solution, 
every  goal  serially.   Go  to  step  4. 


We  have  satisfied 


Step.  4    Form  ISGP  II  -  Tradeoff  Table. 

ISGP  II  -  Tradeoff  Table 


max  f i ( x ) 


m  i  n  f  2  (  x ) 


m  i  n  f  r>  (  x  ) 


f  540.00 

MAG(78S)  441.39 

t~  93.34 

w  1  .00 


8  .44 

30.78 

110.00 

1  .00 


2  .  26 
3.15 
6.  29 
1  .00 


Iteration  1 

Goal  495.33 

PS  412.77 

AS1  495.33 

AS2  412.13 

AS3  412.77 


18  .60 
21.43 
57.22 
18.60 
21  .43 


2  .66 
2.66 

3  .  30 
2.88 
2  .  66 


This  table  is  presented  to  the  DM  for  evaluation. 

§.t£E  5   Evaluation  phase. 

Step  5.1   We  assume  the  DM  is  not  satisfied  with  any  of  the 
solution  sets.  Go  to  Step  5.2. 
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Step_  5.2   Set  new  goals  set,  b. 

We  assume  that  the  DM  realizes  that  unless  all  his/her 
previous  goals  are  relaxed.   no  compromise  solution  will   be 


obtained.   A  new  set  of  goals  is  given  at  80%   of 
i.e.  b  -  (450.67.  28.75.  3.07). 


1 


1 


'I 


Steg  6    Calculate  The  Deviations  Of  The  New  Goals  From  PS. 

t>!  >  fj0  implies  fj"1"  =  37.90  -  450.67  -  412.77 

b2  >  f2°  implies  f2"  =  7.32  =  28.75  -  21.43 

b3  >  f3°  implies  f3~  ■  0.41  -  3.07  -  2.66 

Step_  7   Consistency  check  for  new  b. 

There  are  two  f,"  >  0  when  one  fj+  >  0  and  that  implies  the 
new  goal  set  is  acceptable.  Go  to  step  2.1  and  start  the  second 
iteration  . 

Iteration  No^2 

After  a  new  goal  b  is  given,  the  procedure  in  the  first 
iteration  is  repeated. 
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The  ISGP  II  -  Tradeoff  Table  is  given  as: 


max  f  j (x ) 

540.00 

441 . 39 

93.34 

1  .  00 


min  f 2 (x ) 

8  .44 

30  .78 

110. 00 

1  .00 


min  f  3 ( x ) 
2  .26 
3.  15 
6.29 
1  .00 


MAG(78%) 

f  " 


Iteration  2 
Goal 

PS 
AS1 
AS2 
AS3 


450 .67 
435 .66 
450.67 
435.66 
435 . 66 


28  .  75 
28  .  75 
36.58 
28.  75 
28  .75 


3.  07 
3  .  07 
3.07 
3.07 
3.07 


The    tradeoff    table    is    presented    to    the    DM    for    evaluation. 

Evaluation : 

The  ISGP  II  -  Tradeoff  Table  is  presented  to  the  DM.  The  DM 
finds  that  when  b2  and  b  3  are  satisfied,  bj  receives  less 
attention  than  expected.  Again,  unless  all  goals  are  relaxed, 
there    won't    be    any    agreement    on    the    possible    solutions. 

Set    new    goal s : 

The  new  goal  is  set  at  b  =(  442.00,  30.00.  3.10). 

Calculate  the  deviations  of  the  new  goals  to  PS: 
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bj  >  fj0  implies    f ±+    '    6.36=  442.00  -  435.66 


b2  >  f2°  implies    f2   »  1.25=  30.00  -  28.75 


b3  >  f3°  implies    f3   =  0.03=  3.10  -  3.07 


Consistency  check: 

a  A 

There  are  two    f,~  >  0  when  one    fj  +  >  0  and  that  implies 

the  new  goal  set  is  acceptable.   Go  to  step  2.1  and  start  the 

third  iteration. 

Ii£!iaii°Il  ^0^3 

The  procedure  in  iteration  1  is  repeated.   The  ISGP  II  - 
Tradeoff  Table  is  given  as: 


max  f i ( x ) 

540.00 

441 . 39 

93  .  34 

1  .00 


min  f 2 (x ) 

8.44 

30.78 

110.00 

1  .00 


min  f 3 ( x I 
2.26 
3.  15 
6  .  29 
1  .00 


MAG(78%) 
f" 


w 


Iteration  3 

Goal 

PS 

AS1 

AS2 

AS3 


442 .00 
438 .00 
442 .00 
438.00 
438. 00 


30.00 

30  .00 

31  .67 
30  .  00 
30  .  00 


3  .  10 
3.10 
3  .  10 
3.  10 
3  .  10 
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The  tradeoff  table  is  presented  to  the  DM  for  evaluation. 

Evaluation : 

The  ISGP  II  -  Tradeoff  Table  is  presented  to  the  DM.  Assume 
the  DM  is  satisfied  with  the  PS,  i.e..  the  final  preferred 
solution  is  reached.   Go  to  step  8. 

Steg  8 

The  print-out  of  the  preferred  solution  is  given  as: 

X  -  (3.00.  0,  0.  10.00,  8.00.  4.00) 

f  =  (438 .00,  30.00.  3.10). 
And  the  ISGP  II  algorithm  terminates. 
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4 .4       Remarks: 

1.  The  algorithm  for  ISGP  II  is  an  efficient  procedure  for 
solving  MODM  problems  with  multiple  DMs.  especially  the 
incorporation  of  information  of  f  ,  f  ~ .  MAG,  MAR,  the  principal 
and  auxiliary  solutions.  ISGP  II  can  act  as  a  source  of 
information  for  the  DMs  and  also  help  in  finding  consensus 
agreements  among  DMs  within  a  few  iterations.  The  major  advantage 
of  ISGP  II  is  that  any  single  objective  optimization  technique  can 
be  used  to  obtain  nondominated  solutions.  However,  ISGP  II  can 
only  handle  linear  models  in  which  a  convex  set  is  guaranteed  in 
forming  a  feasible  region.  If  the  feasible  region  is  a  non-convex 
set,  caution  should  be  paid  in  using  ISGP  II.  (See  Fig.  4.2  and 
Fig.  4.3  for  graphical  demonstration  of  convex  and  non-convex 
set.  ) 

2.  If  the  o verachie vement  variables,  d+ ,  are  not  included  in 
expressions  of  ISGP  II.  converting  the  objectives  to  constraints, 
the    formulation    is    given   as 


fj(x)     ♦    dj-     (     fj       -    fj"     ) 


(4.12) 


f,       )    -   bj.       i    €    I 

When  using  this  set  of  formulations  in  the  nutrition  problem,  we 
obtained  an  infeasible  solution.  For  example,  we  intentionally 
selected  a  set  of  b  =  (  94.  110,  6.20  )  and  formed  a  LP  model  as 
f ol lows : 

min    (    £    d_       ) 


124 


Fig.  4.2  Convex  with  two  objectives  in  nonlinear 
problem. 


Feasible 
Region 


Fig.  4.3  Nonconvex  with  two  objectibes  in 
nonl inear  problem. 
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x  e  x 

fjlx)  +  446.66  df  =  94 

f2(x)  -  101.56  d2~  =  110 

f3(x)  -  4.03  d3~  •  6.20 


This  LP  problem  gives  an  infeasible  solution.  In  order  to  amend 
that  pitfall.  d+  should  be  included  into  expression  (4.12)  and 
(4.13)  should  be  formulated  as  follows: 

fj(x)  *    (  if   -  dj+  )  (  fj*  -  fj"  )  =  bj,   j  €  J      (4.13) 

fjlx)  -  (  dj"  -  if    )  (  tf   -  fj*  )  -  bj,   16  1 

The  test  results  for  (4.12)  and  (4.13)  are  given  in  Table  4.3. 

(Notice  that  the  achievement  function   for  both    formulations  is 

i 
min  (   Z  d_   ) ) 

r=i   ' 

Table  4.3   Test  results  for  nutrition  problem 


expression  (4.12) 


expression  (4.13) 


b 

f 

f 

(  94  , 

110  .  6 . 20) 

infeasible 

(396 

52. 

17 

99. 

2 

83) 

(300. 

100 .  5 . 00) 

(300.  108.89.  5.00) 

(396 

52  . 

17 

99  . 

2 

83  ) 

(400. 

20 .  2.50) 

(400 ,  24 .05.  2 . 55) 

(413 

11  , 

22 

97, 

2 

55) 

(500. 

100,  3 .00 ) 

(500 .  100 .  4 . 56) 

(462 

72. 

43 

64  . 

3 

00) 

(539. 

8.50.  2 .30) 

(413  .  11  .  22 .97 .  2.55) 

(413 

11  . 

22 

97  . 

2 

55) 
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It  is  obvious  that  the  results  of  (4.13)  are  at  least  equal  to  or 
better  than  those  of  (4.12).  Cases  such  as  b  =  (  94.  110.  6.20  ), 
b=  (300.  100.  5.00).  and  b  »  (400.  20.  2.50)  confirmed  that 
expression  (4.13)  is  indeed  superior  than  expression  (4.12). 
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Chapter    5.  CONCLUSIONS 

Multiple  Criteria  Decision  Making  problems  are  all  around  us. 
Problems  in  private  and  public  sectors  occur  in  such  a  way  that 
multiple,  conflicting,  noncommensurab 1 e  decision  criteria  are  not 
uncommon.  To  make  decisions  based  on  intuitive  thinking  is  no 
longer  applicable,  hence,  mathematical  programming  techniques  are 
introduced  to  help  the  DM  make  rational  decisions.  The  methods 
that  incorporate  mathematical  programming  techniques  are  also 
known  as  MCDM  solution  methods.  These  methods  are  broadly 
classified  as:  (1)  MADM:  methods  that  deal  with  a  finite  (usually 
small)  number  of  predetermined  alternatives.  (2)  MODM:  methods 
that  deal  with  an  infinite  number  of  alternatives,  and  usually 
requires    design    and    search    efforts. 

The  methods  discussed  in  earlier  chapters  are  TOPSIS-MADM, 
TOPSIS-MODM.  and  ISGP  II.  We  will  discuss  them  briefly  in  this 
section.  (1)  TOPSIS-MADM:  the  idea  of  the  method  is  that  a 
preferred  alternative  should  be  the  one  that  is  closest  to  the 
positive-ideal  point  and  farthest  from  the  negative-ideal  point. 
The  terms  "closeness"  and  "farthest"  are  measured  by  distances. 
The  decision  matrix  must  be  normalized  before  any  comparison  is 
made  because  the  attributes  may  have  different  measurement  units. 
The  parameter,  q,  for  the  normalization  equation  can  be  1.  2.  or 
oo  .  Similarly,  the  parameter,  p.  for  distance  measurement 
function  can  also  be  1.  2.  or  oo.  Different  q  and  p  combinations 
may    result    in    conflicting   solutions.      To    resolve    the    conflict,    the 
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relative  reliability  of  a  distance  measurement,  C,  is  introduced 
to  harmonize  the  final  solutions.  (2)  TOPSIS-MODM:  the  concept  of 
this  method  is  inheritated  from  TOPSIS-MADM.  i.e.  closest  to  the 
p o s i t i v e - i d e a  1  and  farthest  to  the  n  e  g  a  t  i  v  e  -  i  d  e  a  1  .  The 
mathematical  structure  of  this  method  is  adapted  from  the  global 
criteria  method,  i.e.  the  distance-measuring  function.  Although 
the  alternative  that  is  closest  to  the  positive-ideal  as  well  as 
farthest  from  the  negative-ideal  rarely  exists,  we  can  manage  that 
by  introducing  the  importance  index.  H  *  (  H  .  Hn  ),  of  distances 
to    PIS    and   NIS    respectively.    And  Similar    to    that    of    TOPSIS- 

MADM,  the  parameter,  p,  in  distance  function  can  be  1,  2,  and  oo. 
Each  p  value  results  in  a  different  set  of  nondominated  solutions. 
The  DM  may  select  one  or  none  of  the  three  solutions.  (3)  ISGP 
II:  ISGP  II  combines  attractive  features  of  both  Goal  Programming 
and  Interactive  methods.  The  DM  is  actively  involved  in  the 
decision  making  process  while  the  necessary  information  from  the 
DM  is  simplified.  The  DM  gains  insight  into  the  problem  due  to 
the  iterative  nature  of  the  algorithm.  The  cases  with  multiple 
decision      makers      are      handled  by      introducing      the      Maximum 

Achievable  Goals  (MAG)  and  the  Maximum  Achievable  Rate  (MAR).  So 
far,     linear    problems    have    been    successfully   solved   by    ISGP    II. 

The  development  of  MCDM  methods  is  more  and  more  oriented 
twoard  practical  uses.  The  computerized  algorithms  that  allow 
man-machine  "conversation"  have  also  gained  a  lot  of  attention. 
MCDM  methods  are  becoming  more  and  more  powerful  in  their  effort 
to      aid    human    beings    in    making    rational    decisions. 
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Appendix  A. 


Eld  ),  and  the   relative   reliability  rate  C„.  The  program  is  from 


Yoon  [1986]. 


PROGRAM  AA 

DIMENSION  X  (50) 

REAL  NSAM.  NSBM.  NSCM 

WRITE  (5.*)'H0W  MANY  VARIABLE ( FIRST  &  LAST) 

READ  (5.901LS.L  LVC 

WRITE  (5,*) 'HOW  MANY  RANDOM  GENERATION1 
READ  (5.91JPRAN 

WRITE(5.200) 
WRITEI5  ,210)1,3,  LVC 
WRITE(5 .220IPRAN 
WRITEO . 200 
WRITEO .2101LS.LVC 
WRITEO.  220)PRAN 
WRITE(3,100) 
WRITE(5 . 100) 
WRITEO,  110) 
WRITEO. 110) 

DO  10  LV  =  LS.LVC.l 
AK=0 .0 
BK  =  0  .  0 

A-l 


CK=0 .0 
SA=0.0 
SP  =  0. 0 
SC=0.0 
SAM  =  0  .0 
SBM=0 .0 
SSM=0 .0 
ZK=0.0 
QK  =  0 .0 
PM=0.0 
NSAM=0. 0 
NSBM=0. 0 
NSCM=0.0 

DO  15  ZK-1 .0 .PRAN. 1 .0 
0K=QK+1 . 0 
AK=0.0 
BK  =  0 .0 
CK=0. 0 
PB=0.0 

DO  20  H=l,LV.l 
X(M)=RAN(I , J) 
20       CONTINUE 
M  =  l 

P-X(M) 

DO  30  M=2 ,LV, 1 
IF(P.LT.XIM)  )  GO  TO  35 
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p=p 

30  ID  30 
35       P  =  X  (M) 
30        CONTINUE 

DO  40  M=l ,LV. 1 

AK=AK+X(M) 

PB-PB+X(M)**2 
40        CONTINUE 
DUAO:  [USERIE.KY175]AA.F0R:59 

BK=PB**. 5 

CK  =  P 

SA=SA+AK 

SB=SB+BK 

SC=SC+CK 
15        CONTINUE 

SAM=SA/QK 

SBM=SB/QK 

SCM=SC/QK 

PM=SAM+SBM+SCM 

NSAM=SAM/PM 

NSBM=SBM/PM 

NSCM=SCM/PM 

WRITE(3,400)LV.SA.SB.SC,SAM.SBM,SCM,NSAM.NSBM.NSCM 

WRITE(5.400)LV.SA.SB.SC.SAM,SBM.SCM.NSAM.NSBM.NSCM 
10        CONTINUE 

WRITEI3, 110) 
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WRITE(5, 1101 
91        F0RMAT(F12.  1  ) 
90       F0RMAT(2IB) 
100      F0RMAK///5X,  '#  OF  VAL  '  .  9X  ,  '  SA  '  ,  14X  ,  '  SB  '  ,  14X  .  '  SC  '  ,  1 IX  . 

S   'SAM' .7X. 'SBM' .7X. 'SCM.7X. ' NSAM ' .6X. 'NSBM' ,6X, 'NSCM' ) 
110      FORMATOX,  120(  '-'  )  ) 

200      F0RMAT(//12X, ' *  SIMULATION  CONDITION  *') 
210      FORMAT (8X,  'NUMBER  OF  VALIABLES  -  FROM '  . 1 3 .  ' TO '  .  13) 
220      FORMAT (3XM' NUMBER  OF  RANDOM  GENERATION  TIME  -'.  F3.0) 
4  00      FORMAT (/5X. I3.GX.3F15.4.3X.GF10.4) 

STOP 

SND 
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Appendix  B 

THE  METHOD  OF  GLOBAL  CRITERION 
The  method  of  global  criterion  Is  to  obtain  an  optimal  vector  that 
minimizes  some  global  criterion:  for  example,  the  sum  of  the  squares  of  the 
relative  deviations  of  the  criteria  (objective  functions)  from  the  ideal 
points. 

The  method  is  used  to  obtain  a  solution  for  the  MODH  problem  formulated 
below: 

max/min  [fjlx).  f2<x).  fglx) fk(-)]  (311 

s.t.   x  6  X   where  I  •  (t  |  |.  (x)  {*  •  S}  0,  1  =  1.2 m) 

where     f  (x):   benefit  objective  for  maximization,  j  6  J 
f  (x) :   cost  objective  for  minimization.  1  €  I 

The  method  of  global  criterion  for  the  vector  problem  of  (3.1)  solves  the 
problem  given  by  (3.2): 

,r  ft  -  'i  w   ]  p  I"  fi(-xl  -  fI  ]  P 

min    F  =  Z       w.   --■!, i +  =   »,   ~" i"    <3-2> 

(X)     P   J6J   J  [   f j  -  f J     j     i€X     I    fT  -  f,  \ 

s.t.    x  €  X. 
p  -  1.2.a 
where  w  .  i  =  1.2 k  are  the  relative  importance  (weights)  of  objectives. 

Boychuk  and  Ovchinnikov  [1973]  have  suggested  p  =  1.  and  Salukvadze 
[1972.1974]  has  suggested  p  =  2  for  the  global  criterion.   Although  global 
criteria  p  can  be  any  value  between  1  and  a,    distances  with  p  =  1,  2.  and  oo 
are  especially  operationally  Important.   Distance  d  (the  Manhattan  distance. 

or  city  block  distance)  and  d  (the  Euclidean  distance)  are  the  longest  and 

shortest  distances  in  the  geometric  sense:  distance  d   (the  Tchebycheff 

oo 

distance)  is  the  shortest  distance  In  the  numerical  sense.  The  'best' 
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solution  chosen  for  the  MODM  will  differ  greatly  according  to  the  criteria 
chosen.   These  criteria  locate  points  which  are  geometrically  different  from 
one  another.   The  problem  with  this  approach  is  to  determine  which  p  would 
result  in  a  solution  that  is  the  most  satisfactory  solution  to  the  Decision 
Maker  (DM).   It  is  also  possible  that  whatever  p  is  chosen,  (3.2)  might  give  a 
solution  which  has  a  particular  objective  value,  f^fx),  that  is  totally 

unacceptable  to  the  DM. 

Case  1.   p  ■  1:  f(x)  and  g(x)  are  linear  functions. 

When  all  the  objectives.  f,(x),  i  =  1.  2 k,  and  all  the  constraint 

functions,  g. (x) .  1  =  1 m.  are  linear  functions,  and  p  =  1.  then  (3.2) 

becomes  a  Linear  Programming  (LP)  problem.   The  simplex  method  for  LP  can 

solve  the  problem. 

Case  2.   p  =  2:  f_(x)  and  jj(x)  are  linear  functions. 

When  both  f(x)  and  g(x)  are  linear  functions  and  p  =  2.  (3.2)  becomes  a 
convex  programming  problem;  more  precisely,  a  quadratic  programming  problem. 

As  is  well  known,  in  a  convex  programming  problem,  a  local  minimum  is 
absolute  (or  global).   The  convex  programming  problem  can  be  solved  by  using 
any  method  of  finding  a  local  extreme,  including  the  gradient  method, 
quadratic  programming,  or  nonlinear  programming  methods. 

The  case  where  p  -  3.  4.  ...etc  are  normally  quite  similar  to  the  case 
where  p  =  2.   We  will  now  consider  the  case  where  p  -  »  and  show  how  it  can  be 
reduced  to  an  LP  problem. 
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Case  3.   p  =  o:  f_(x)    and  £(x)  are  linear  functions. 

r*       io>        r        *ioo 

mill   F   =  Z       w,   !— ig i~~—    '    +  E   w.  j  —  --- r— -  \ 

"   iej  J  !  f  -  f"   '    lei   *  '   f"  -  f 

JfcJ       t   rj    'j     j      1C1      L     i     1     J 
s.t.     X  6  X  (3.3) 

(3.3)  can  be  transformed  to  (3.4)  for  solution  of  the  LP  problem  (see 
Note  (4)  for  this  transformation): 

min   d  (minimizing  deviation) 

00 

s.t.    x  6  X 

I   f1  "  f1,s) 
»4    —  » _ s  d      (^  j  6  J)  (3.4) 

J   |       ff     _  f  '  ■ 

L     J     j     J 

r        *   i 

I  fi(*'  "  fi 
».   -" «-—-   s  <»     (v  1  6  I) 

L   ri    ri      J 
This  is  a  MINIMAX  approach  problem.   Among  all  d  terms  in  the  constraints  we 

are  to  minimize  the  maximum  term  of  d  . 

Case  4.  £(x)  and  g(x)   are  nonlinear  functions. 


If  any  one  or  all  of  the  objectives.  f.(x).  1  =  1.2 k.  and  the 

constraints.  g,(x).  i  »  1.  2 m.  are  nonlinear  functions.  (3.2)  is  a 

nonlinear  programming  problem  (for  both  p  •  1  and  p  =  2).   For  p  =  <o  and  one 
or  more  than  one  of  f,(x)  and/or  g.(x)  are  nonlinear  functions.  (3.4)  format 

solves  the  problem,  but  it  is  a  nonlinear  programming  problem.   A  variety  of 

nonlinear  programming  methods  (e.g.  the  sequential  unconstrained  minimization 

technique  (SUMT)  [Hwang,  et  al . .  1975],  or  the  generalized  reduced  gradient 

method  (GRG)  [Abadie  &  Carpenter.  1969;  Hwang  et  al . .  1972])  can  be  used  for 

solving  the  problem.   A  new  approach  which  is  called  iterative  nonlinear  goal 
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programming  is  an  effective  method  for  solving  this  problem  (Hwang  S  Masud. 
1979:  see  Sec.  2.2.2.D). 

NUMERICAL  EXAMPLE  1.   (Production  Scheduling  of  the  Hardee  Company) 

The  Hardee  company  makes  two  kinds  of  dolls.   Doll  A  is  of  higher  quality 
and  Doll  B  is  of  lower  quality.   The  respective  profits  are  $0.40  and  $0.30 
per  doll.   Each  Doll  A  requires  twice  as  much  labor  as  a  Doll  B,  and  if  all 
dolls  were  of  type  B.  the  company  could  make  500  per  day.   The  supply  of 
material  is  sufficient  for  only  400  dolls  per  day  (both  A  and  B  combined). 
The  problem  assumes  that  all  the  type  A  and  type  B  dolls  the  factory  can  make 
could  be  sold,  and  that  the  best  customer  of  the  company  wishes  to  have  as 
many  as  possible  of  type  A  doll.   The  manager  realizes  that  two  objectives: 
(1)  the  maximization  of  profit,  and  (2)  the  maximization  of  production  of  Doll 
A,  should  be  considered  in  scheduling  their  production. 
This  VMP  can  be  formulated  as  follows: 
max  f,(x)  •  0.4x  +  0.3x2 


subject  to 


max  f„(x)  =  x 


gjfx)  =  Xj  +  xg  £  400  "I 


g2(x)  -  2Xj+  x2  £  500  >  x  €  X 

where  x  and  x  are  the  number  of  Doll  A  and  Doll  B  produced  respectively. 

The  steps  of  the  method  of  global  criterion  are  as  follows: 
Step  0.   Construct  the  Positive-Ideal  Solution  (PIS)  and  Negative-Ideal 
Solution  (NIS)  payoff  tables. 

It  involves  solving  the  objectives  as  individual  LP  problems: 
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0.1  (PIS)      1)    max  fjfx)  =  0.4XJ  +  0.3x2 

s.t.    x  G  X 
The  solution  Is:   Xj  =  100,  xg  =  300;   fj(x*)  =  130 

Point  B  is  the  ideal  solution  for  f.(x)  (See  Figs.  3.1  and  3.2). 

Similarly  for  the  second  objective: 
2)    max  f2(x)  -  xx 

s.t.  x  6  X 


The  solution  is:  Xj  =  250,  x2  =  0;   fgfx*)  ■  250 
Point  C  is  the  ideal  solution  for  f,(x>- 
The  PIS  payoff  table  is  therefore: 


fl      f2 


f   I  130*     100   I  100    300 
f   I  100-     250*  I  250      0 


Note  that  in  this  payoff  table,  row  1  gives  the  solution  vector  x*  =  (100, 
300)  which  maximizes  the  objective  f^x*)    •   130:  fg  (100,300)  =  100  is  the 

value  taken  on  by  the  objective  f  when  f  reaches  its  maximum.  Similarly. 

row  2  gives  the  solution  vector  x*  =  (250,  0)  which  maximizes  l^   (x*)  =  250; 

f,  (250,  0)  =  100  is  the  value  taken  on  by  f  when  f  reaches  its  maximum. 
1  18 

The  PIS  is  point  P.  (f  .  £  )  =  (130,250)  shown  in  Figure  3.2.   This  is  a 

Utopia  solution,  which  is  not  feasible. 

0-2.   (NIS)    1)     min  f^x)    -  .4x1   +  .3x2 

s.t.    x  €  X 
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The  solution  Is:       Xj  =  0.   xg  =  0:   fjlx-)  "  ° 
2)    rain  f2(x)  =  *1 
s.t.   x  €  X 


The  solution  is:    Xj  =  0,  0  S  «2  i  400:    f2(--)  =  ° 


The  NIS  payoff  table  Is  therefore: 


fl      f2   I   Xl 

I 


fl|   °"      ° 

f  I  120      0-0     0  i  x2  S  400 

The  NIS  is  point  N.  (f~.  f~_)  =  (0,0).  shown  in  Figure  3.2. 


Step  1.  Obtain  the  preferred  solution. 

A  preferred  solution  is  a  nondominated  solution  which  is  a  point  on  the 
segment  of  straight  line  BC.   (In  Figs.  3.1  and  3.2  all  points  on  the  segment 
of  straight  line  BC  are  nondominated  solutions).  Assume  equal  weights  for  the 
two  objectives  (w  ■  w  -  0.5). 

Case  1.   p  =  1 

The  preferred  solution  for  p  =  1  in  (3.2)  is  to  find  x  =  (x, .»_)  so  as  to 

r  130  -  (0.4x  +  0.3x  )   1       r  250  -I,  1 
rain  Fp=1  =  0.5*  I  jjg-^ |  ♦  0.5*1  -555—5-  I 

=  1  -  0.00354X   -  0.00115X 

S.t.    X  €  X 
This  is  a  linear  programming  problem.   The  solution  is  given  by 
x  =  [250,  0] 

f  •  [100,  250] 
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which  is  point  C  in  Figs.  3.1  and  3.2. 
Case  2.   p  =  2 

The  preferred  solution  for  p  =  2  in  (3.2)  is  to  find  x  =  (x,,  xg)  so  as 

to 

2  2 

r   130  -  (0.4x   +  0.3x  )   1  f  250  -  x   1 

•in  Fp=2  -  0.5*  I  jjjJ -2-  J   ♦  0.5*  I  — -2B5— -  J 

s.t.   x  e  x 

This  is  a  nonlinear  programming  problem.   The  sequential  unconstrained 

minimization  technique  (SUMT)  [Hwang,  et  al .  1975]  is  used  for  solving  the 

problem.   The  solution  is 

X  =  [230.70.  38.60] 

f  ■  [103.90.  230.70] 

which  is  point  D  in  Figs.  3.1  and  3.2. 

Case  3.   p  ;  oo 

The  preferred  solution  for  p  =  oo  in  (3.2)  is  to  find  x  so  as  to 

min   d 
oo 

s.t.   x  €  X 

r  130  -  (.4x  +  .3x  )  T 
0.5    ,  „„i.„. ..?.  |  s  ^ 

r  250  -  x   1 
°-5    I  -250--0--  J   *  d» 

This  is  a  LP  problem.   We  can  use  the  LP  simplex  method  to  solve  it.   The 
solution  is: 

x  =  [210.  84] 

f  =  [109.2.  210] 
which  is  point  E  in  Figs  3.1  and  3.2. 
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As  shown  in  Fig.  3.1  for  p  =  1.  the  function  to  be  minimized.  F   .  is  a 

linear  function:  therefore  the  preferred  solution  is  determined  by  the  slope 
of  line  BC  and  that  of  F  ,  (i.e..  min  F  .  =  1  -  0.00354  x  -  0.00115  x  is 

p=l  p=l  1  c, 

equivalent  to  max  (-F   )  =  0.00354  x  +  0.00115  x„  -  1.)  Since  the  slope  of 

F  .  is  less  than  that  of  line  BC  (the  line  representing  the  nondominated 
p=l 

solutions)  ,  point  C  is  the  preferred  solution.   For  p  =  2.  the  function  to  be 

minimized,  F   „.  is  a  nonlinear  (quadratic)  function,  and  the  minimum  point 
p=2 

for  F  „  is  point  D  in  Fig.  3.1.  which  is  the  point  where  the  F    curve 
p=2  P-c 

touches  the  segment  of  line  BC.   For  p  ■  ».  the  function  to  be  minimized, 

F   ,  cannot  be  drawn  on  Fig.  3.1  but  can  be  represented  in  Fig.  3.2  as  point 
p=oo 

E.   Note  that  point  E  is  the  point  that  the  largest  square  drawn  between  PIS 

and  BC  touches  line  BC  in  a  normalized  objective  function  space 

representation. 

The  DM's  preferred  solution  for  this  example  should  be  one  of  the 

nondominated  solutions  (represented  by  points  on  the  segment  of  straight  line 

between  B  and  C  in  Figs.  3.1  and  3.2). 

The  solutions  for  F   ,.  point  C,  and  for  F  ,.  point  D.  and  F   .  point 
p=l  r  p=2  p-<» 

E,  derived  by  this  method  may  not  satisfy  the  DM  at  all. 

NUMERICAL  EXAMPLE  2.   (The  Nutrition  Problem) 

The  nutrition  problem  is  to  find  the  quantities  of  certain  foods  that 
should  be  eaten  to  meet  certain  nutritional  requirements  at  a  minimum  cost. 
Information  on  six  foods  is  given  in  Table  3.1. 

The  problem  is  to  find  the  quantities  of  these  foods  that  satisfy:   (i) 

the  minimum  daily  requirement  of  vitamin  A  and  iron:  and  the  balanced  daily 

requirement  of  food  energy  and  protein;  (ii)  the  minimum  cholesterol  daily 

intake;  (ill)  the  maximum  carbohydrate  daily  intake:  and  (iv)  the  minimum 

cost. 
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The  maximum  daily  intakes  of  the  foods  are:  6  pints  of  milk,  1  pound  of 
beef,  1/4  dozen  eggs.  10  ounces  of  bread,  10  ounces  of  lettuce  and  salad,  and 
4  pints  of  orange  Juice. 

Let  the  daily  diet  requirements  of  milk  (in  pints),  beef  (in  pounds), 
eggs  (in  dozens),  bread  (in  ounces),  lettuce  and  salad  (in  ounces),  and  orange 
juice  (in  pints)  be  x  .  x  ,  x  ,  x  .  xg,  and  xg .  respectively.  Then,  the  HODM 

problem  is  to  find  x  =  (x  .  x Xg)  so  as  to: 

max   f  (x)  =  24x  +  27x  +  15x4  +  1 . lxg  +  52xg  (carbohydrate  intake) 

min   f„(x)  =  10x.  +  20x„  +  120x„  (cholesterol  intake) 

min   f3(x)  =  0.225  Xj  +  2.2x2  +  0.8x3  +  0.1x4  +  0.05x5  +  0.26Xg(cost) 

s.t.   g  (x)  -  720x  +  107x2  +  7080x3  +  134x5  +  lOOOXg  £  5000 

g  (x)  =  0.2XJ  +  lO-lXg  +  13-2x3  +  0.75x4  ♦  0.15x5  +  l-2Xgi  12.5 

g  (x)  =  334x  +  1460X  +  1040x3  +  75x4  +  17.4xg  +  240xgi  2500 

g  (x)  =  18XJ  +  151x2  +  78x3  +  2.5x4  +  0.2xg  +  4xg  £  63 


g5(x)  =  Xj  S  6.0 

g6(x)  =x2  i  1.0 

g?(x)  =  x3  S  0.25 

88(5)  -  x4  *  100 

gg(x)  =  xg  £  10.0 

g1Q(x)  =  x6  U.O 


x  i  0,   i=1.2. 


Let  the  constraints  be  denoted  by  x  €  X. 
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Table  3.1  Nutrition  and  cost  Information  of  six  foods 


Recommended 
Milk   Beef   Eggs   Bread  Lettuce  Orange      Daily 

salad    juice    Allowance 
(pint)  (pound) (dozen) (ounce)  (ounce)   (pint)    for  adults 


Vitamin  A  (I.U. 

Food  energy 
(calories) 

Cholesterol 
(unit) 

Protein  (g) 


720    107 
344   1460 


10 


18 


20 


151 

Carbohydrate  (g)   24     27 
Iron  (mg)  0.2   10.1 

Cost  ($)  0.225  2.2 


7080     0       134     1000      5000 
1040    75       17.4    240      2500 

120     0        0        0 


78           2.5 

0.2 

4 

0         15 

1.1 

52 

13.2      0.75 

0.15 

1.2 

0.8     0.1 

0.05 

0.26 

12.5 


Step  0  Obtain  the  ideal  solutions  and  construct  payoff  tables. 

Using  the  same  procedure  as  in  Example  1.  we  obtain  all  the  PIS's  and  NIS's. 

The  payoff  tables  are  given  as  follows: 

Positive  Ideal  Payoff  Table 

MAX  Carb  MIN  Choi  MIN  Cost      Milk   Beef   Eggs  Bread  Lettuce   0J 
MAX  Carb  540. 00«  80.00    6.09      6.00   1.00   0.00   10.00   10.00   4.00 
MIN  Choi   380.39    8.44*    3.47      0.00    0.42    0.00   10.00   10.00   4.00 

MIN  Cost   281.67    67.83     2.26*     3.78  0.00    0.25   10.00    0.00    0.79 


Negative  Ideal  Payoff  Table 

MAX  Carb  MIN  Choi  MIN  Cost 
MAX  Carb   93.34"  71.53    3.08 
MIN  Choi   321.00   110.00"   4.75 
MIN  Cost   540.00   110.00     6.29" 


Milk  Beef  Eggs  Bread  Lettuce  0J 

2.50  0.83  0.25  0.00  10.00   0.00 

6.00  1.00  0.25  10.00  10.00   0.00 

6.00  1.00  0.25  10.00  10.00  4.00 
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Step  1  Obtain  the  preferred  solution.   (For  w  =  w  ■  w  =  0.33) 
Case  1  p  ;  1 : 


min   F   .  =  0.33 
p-1 


0.33 


540  -  (24x  +  27x  +  Ox  +  lSx  +  1 . lx5  +  52xg  )1 
__________ 


(10X   +  20x2  +  120x3)  -  8.44   1 
---------- 


(.225x  +  2.2x2  +  .8x3  +  . lx4  +  .05x.  +  . 26xg )  -  2.26  T 
___________  j 

►  .6824x„  +  1.3801X.  -  .0088x„  +  ,0099x_ 


-  ,0519x.  +  .5651 

D 

s.t.   x  6  X 
This  is  a  LP  problem.   The  solution  is 

X  =  [3.76.  0.01.  0.13.  10.00.  10.00.  0] 
f  =  [253.89.  55.51.  2.65] 
Case  2  p=2: 


min  F   „  =  0.33 
p=2 


r  540  -  (24x   +  27x   +  Ox   +  15x   +  1 .  lx  +  52x  )   1 
.  1      £•  J      y       o     o 


r   (10x  +  20x  +  120x  )  -  8.44  1 
+  0.33   |  — ----- ---?-------3 , 


r   (.225x.  +  2.2x„  +  .8x„  +  .lx,  +  .05xr  +  .26x„)  -  2.26   1 
+  0.33  |  ---i---?--?------^------§_--e , 

s.t.     x  €  X 

This  is  a  nonlinear  programming  problem.  The  Sequential  Unconstrained  Minimi- 
zation Technique  (SUMT)  is  used  to  solve  this  problem.   The  solution  is 
x=  [2.48,  0.13.  0.00.  9.94.  1.08.  4.00] 

f  =  [420.93.  27.75.  2.92] 
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Case  3.   p  ;  a 


min    d 

9 


s.t.   x  6  X 

T  540  -  (24x  +  27x   +  OXg  +  15x4  +  1 . lXg  +  52Xg)  1 
0.33  |  540-"93734 ""  \      *  d. 

T   (10x,  +  20x   +  120x  )  -  8.44   1 

o.33 1  — i — ni---8T4i5 J  *  d» 

r   (.225x   +  2.2x   +  .8x3  +  . lx4  +  .05Xg  +  -26x6)  -  2.26  1 

0.33      |  6729"~2726  S     oo 

L  J 

This   is  a  LP  problem.      The  solution  is: 

X  =    (3.08.    0,    0.    10.00,    8.63,    4.00) 

f  -    (441.43.    30.81.    3.15) 

The  solutions  of  the  three  cases  are  presented  to  the  DM.   The  DM  may 

pick  any  one  of  them,  or  may  not  agree  with  any  of  them. 

Notes: 

(1)  The  solutions  of  the  nutrition  problem  are  summarized  in  the  following 

table.   Also  in  the  table  is  a  set  of  solutions  for  a  given  set  of 

sights  [0.30.  0.20.  0.50]  for  [fjlx).  fjfx).  f 3 (x) ] .   It  is  obvious 


wei 


that  when  the  weight  for  f,(x).  the  cost  objective,  is  0.50  rather  than 

0.33.  the  cost  is  reduced  to  2.32  from  2.65  for  p=l;  to  2.65  from  2.92 
for  p=2,  and  to  2.8  from  3.15  for  p=o. 
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f! 

f2 

f3 

Xl 

X2 

X3 

X4 

X5 

X6 

Positi 

ve  ideal  solution 

540.00 

* 
8.44 

* 
2.26 

Negative  ideal  solution 

93.34" 

110. 00~ 

6.29~ 

(0.33.0.33.0.33) 

253.89 

55.51 

2.65 

3.76 

0.10 

0.13 

10.00 

10.00 

0 

P=l 

(0.30.0.20.0.50) 

249.91 

70.85 

2.32 

3.98 

0.05 

0.25 

10.00 

2.65 

0 

(0.33.0.33.0.33) 

420.93 

27.75 

2.92 

2.48 

0.13 

0 

9.94 

1.08 

4.0C 

0=2 

(0.30.0.20.0.50) 

417.70 

29.15 

2.65 

2.48 

0 

0.03 

io.oo 

0.60 

3.9S 

(0.33.0.33.0.33) 

441.43 

30.81 

3.15 

3.08 

0 

0 

10.00 

8.63 

4.0C 

p=oo 

(0.30.0.20.0.50) 

440.48 

34.37 

2.80 

3.44 

0 

0 

10.00 

0 

4.0C 

(2)  There  can  be  infinitive  solution  sets  for  a  M0DM  problem  solved  by 

the  method  of  global  criteria,  because  the  global  criteria  p  can  be  any 

integer  between  1  and  <».  Among  the  infinitive  p's.  Yoon  [1986]  mentioned 

three  meaningful  p's.  i.e.:   p  =  1.  2.  and  oo.   For  p  ■  1.  distance  between  two 

points  is  the  longest  distance  in  the  geometric  sense.   For  p  ■=  2.  distance 

between  two  points  is  the  shortest  distance  in  the  geometric  sense.   For  p  = 

o.  distance  between  two  points  is  the  shortest  distance  in  the  numerical 

sense.   Three  parameters,   p  ■  1.  2.  and  oo  are  widely  used  to  solve  the  M0DM 

problem. 

Since  the  solutions  will  be  different  among  different  p's.  we'll  face 

three  different  solutions.   To  choose  a  preferred  solution  depending  on  human 

cognitive  power  is  not  scientific  nor  reliable.  Yoon  [1986]  introduces  the 

relative  reliability  of  a  distance  (d  )  measurement  to  reconcile  several 

P 

compromise  solutions. 

The  relative  reliability  of  d  is  written  as: 

P 

X  -  Eld  )/{E(d,)  +  Eld  )  +  E(d  )} 
p     p      1       £  00 
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where  each  X  =  E(d  )'s,  (mean  of  d  1,  Is  obtained  using  the  Monte  Carlo 
P      P  P 

simulation  technique. 

The  three  distance  measurements  (d  )  are  combined  by  letting  Xp  be 

weights  in  a  linear  combination  form.   The  unified  distance  measurement  (d) 
for  each  individual  alternative  is  calculated  by: 

d  "  \  dl  +  X2  d2  +  X»  d« 
Each  alternative  has  a  d  value.   We  compare  the  d's  and  choose  the  least  one 
as  our  best  solution. 

(3)  Advantages  and  disadvantages  of  this  method. 

The  advantage  of  this  method  is  that  in  the  process  of  obtaining  the 
solution,  the  DM  will  not  be  disturbed  by  the  analyst,  which  is  preferred  by 
the  DM. 

The  disadvantage  is  that  the  analyst  makes  all  necessary  assumptions 
about  the  DM's  preference.   This  is  a  difficult  task  even  for  well  trained  and 
knowledgeable  analysts. 

(4)  We  will  use  a  numerical  example  to  illustrate  that  the  solution  to 
equation  (3.4)  is  equivalent  to  the  solution  for  p  =  »  for  the  Method  of 
Global  Criterion. 

The  PIS  and  NIS  of  the  nutrition  problem,  and  an  arbitrary  set  of 
nondominated  solutions  is  given  as: 

Numerical  Example  2.  (Nutrition  Problem) 

Max  Carb  Min  Choi        Min  Cost 
PIS      540.00  8.44  2.26 

NIS 93.34 110.00 6.29 

f     412.55  20.43  2.75 

Milk   Beef    Eggs    Bread   Lettuce    OJ 
X  =  [   2.04   0.00     0.00     10.00      5.00    4.00] 
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If  w     =  w     =  w      is   assumed: 


„  f     540   -  412.55      ]    *      \     20.43   -   8l44_  ]    *      [  JiZLZ-SiS-  ]    " 

mln  Foo  =    '    -Ko"=~557§4"  j  [  "TloTo   -  8.44      \  [     6.29   -   2.26     J 


CO 


=  [0.2853]™  +  [0.1181]™  +  [0.1216] 
For  p  =  50  or  p  ■  100  to  replace  p  -  co,  we  set: 

—9ft  —47  —46 

F      =  5.822  x  10     *  4.026  X  10     +  1.756  X  10 
p=50 

F      -  3.390  X  10"55  ♦  1.679  X  10"93  ♦  3.114  X  10~92 
p=100 

For  both  F     and  F  .„„.  from  the  order  of  magnitude  of  point  of  view. 
p=50      p=100 

the  first  term  of  the  summation  dominates,  and  the  second  and  third  terms  can 

be  neglected.   Therefore,  the  p  =  to  problem  is  equivalent  to  minimization  of 

the  maximum  dominated  term.  i.e..  MINIHAX  approach. 
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ABSTRACT 
Decision  making  problems  which  involve  multiple,  conflicting, 
and  noncommensurable  criteria  are  frequently  encountered  in  both 
private  and  public  organizations.  Currently,  decision  problems 
are  solved  by  managements  using  their  intuitive  judgements.  Such 
a  decision  process  is  dangerous  and  inconsistent.  To  help  the  DM 
make  reliable  decisions,  researchers  have  poured  their  efforts 
into  forming  various  solution  methods  for  all  kinds  of  MCDM 
problems.  The  MCDM  problems  can  be  broadly  categorized  into  (1) 
multi-attribute  problems,  and  (2)  multi-objective  problems.  Each 
category  has  different  problem  characteristics.  None  of  the 
solution   methods    are    capable    of    solving   both    kinds    of    problems. 

In  this  study,  this  author  discussed  three  solution  methods, 
namely,  TOPSIS-MADM,  TOPS  I S - MODM ,  and  ISGP  II.  The  first  one 
deals  with  the  mu 1 t i -a t t r i bute  problem.  It  utilizes  a  distance 
measurement  function  and  attribute  normalization  function  to 
obtain  an  preferred  solution  which  is  closest  to  the  ideal 
solution  and  farthest  to  the  negative  ideal  solution.  The  second 
one  solves  mu 1 t i -ob j ec t i ve  problems.  It  is  an  extension  of  the 
method  of  Global  Criteria  and  has  properties  inheritated  from 
TOPSIS-MADM.  The  relative  importance  factor.  H,  is  introduced  to 
combine  two  distance  measurements:  distance  to  positive  ideal 
solution,  and  distance  to  negative  ideal  solution.  The  last  one 
solves  multi-objective  problems  too.  It  combines  attractive 
features  of  the  STEM  method,  the  concept  of  Goal  Programming,  and 
an    iterative    approach    to    systematically    derive    a    set    of    preferred 


I 


solutions.   The  DM  is  actively  involved  in  the  solution-finding 
process . 

The  formations  of  MCDM  solution  methods  provide  the  DMs  with 
systematic  approaches  to  derive  better  decisions.  Computerization 
of  those  solution  methods  enhances  the  power  of  the  methods  while 
reducing  the  burden  of  calculations,  hence  makes  MCDM  solution 
methods  more  realistic  in  real-world  applications. 


